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ABSTRACT 


Measurements  of  the  ratio  of  the  neutron-proton 

electric  form  factors,  G„  /G„  ,  were  made  from  elastic 

'   E    E  ' 
n   p 

electron-deuteron  scattering  to  a  precision  of  approx- 

2 
imately  1  to  2%    for  the  range  of  momentum  transfers,  q  , 

2         -2 

of  0,10  <  q   <  0,8  F   ,  and  for  electron  scattering  angles 

of  ^5°  to  120°.   It  was  found  that  within  experimental 

errors  the  slope  as  determined  from  the  ratio  Gp  /GF  , 

n   p 
agrees  with  the  thermal  neutron-electron  interaction  slope 

when  relativistic  corrections  and  proper  deuteron  wave 

functions  are  applied  to  the  electron-deuteron  scattering 

results . 
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I.   INTRODUCTION 


A.   EXPERIMENTAL  OBJECTIVES 

The  charge  form  factor  of  the  neutron,  G„  ,  as  deter- 

n 
mined  by  electron  scattering  is  found  to  be  essentially 

2  2 

zero  in  the  range  of  momentum  transfers,  q  ,  of  0.3  <_  q 

_2 

<_  100  F   .   This  fact  is  in  apparent  conflict  with  the 

results  obtained  by  the  scattering  of  thermal  neutrons 
from  atomic  electrons.   By  methods  to  be  described  later, 

measurements  of  the  slope  of  the  neutron  charge  form  factor 

2      -9   -2 
at  very  low  momentum  transfers  of  q   -  10   F   ,  show  that 


dGn 


n 


d(q2) 


Jq2+0 


=    +0.0193    ±    0.0004    F' 


(1-1) 


It  is  very  doubtful  that  a  non-zero  slope  for  GF   at  q  ->0 

2        -2   n 
can  be  reconciled  with  GF   -  0  at  q   >  0.3  F 

n 

Drickey  and  Hand  [1]  in  1962  measured  Gp  /G„   for 

2  2-2  n    P 

values  of  q  ,  0.3  <  q   £  2.2  F   to  an  accuracy  of  about  2%, 

Their  results  yielded  G^   -  0,  but  when  reanalyzed  by  Casper 

n 
and  Gross  [2]  (1966)  using  recently  developed  deuteron 

relativistic  corrections  and  new  deuteron  wave  functions, 

the  value  of  Gp   was  brought  into  fair  agreement  with  the 

n 
value  of  (1-1).   Thus  the  entire  problem  is  further  compli- 
cated by  the  extreme  sensitivity  of  GF   to  the  different 

n 
deuteron  models  assumed.   Hence,  if  very  accurate  measure- 

2 
ments  of  GF  ,  for  q   small,  can  be  made  by  electron-scattering 

n 
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techniques,  a  comparison  with  the  value  (1-1)  could  be  made, 

which  might  yield  an  answer  to  the  "correct"  deuteron  model 

to  be  chosen. 

It  is  the  purpose  of  this  experiment  to  (a)  demonstrate 

the  basic  agreement  between  the  thermal  neutron-electron 

slope  value  of  GE   and  the  electron  scattering  G„  ;  and 

n  Jn 

(b)  try  to  determine  which  deuteron  wave  function  (model), 

best  describes  the  deuteron. 

B.   NUCLEAR  STRUCTURE  PROM  ELECTRON  SCATTERING 

Since  the  early  1950 's  production  of  beams  of  high 
energy  electrons  by  linear  accelerators  (LINACS)  has  made 
possible  the  study  of  the  electromagnetic  structure  of  the 
nucleon.   It  is  expected  that  the  dimensions  of  any  particle 
of  mass  m  is  of  the  order  of  h/mc ,  which  for  a  nucleon  is 
about  2  x  10     cm.   Hence,  to  investigate  the  structure  in 
this  size  range  it  is  necessary  to  use  a  probe  having  about 
these  wave  lengths  or  shorter.   Linacs  produce  electrons 
with  just  such  wave  lengths.   The  100  MeV  NPS  Linac  thus 
"sees"  only  the  outer  structure  or  surface  of  the  nucleus. 

Electron  scattering  has  several  definite  advantages 
over  other  methods  of  nuclear  structure  research.   Briefly, 
the  important  reasons  that  electron  scattering  is  such  a 
useful  tool  are  the  following: 

(a)   The  basic  interaction  between  the  electron  and 
target  nucleus  is  known;  the  electron  interacts  with  the 
charge-current  density  and  magnetization  density  of  the 
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nucleus.   Since  the  electromagnetic  interaction  is  rela- 
tively weak,  of  the  order  of  a  =  1/137,  the  fine  structure 
constant,  measurements  on  the  target  nucleus  can  be  made 
without  greatly  disturbing  its  structure. 

This  is  in  contrast  to  the  situation  with  protons, 
neutrons  and  alpha  particles,  since  they  are  strongly  inter- 
acting particles,  and  the  scattering  mechanism  cannot  be 
clearly  separated  from  structure  effects  in  the  target. 

(b)   Photons  interact  electromagnetically  with  the 
nucleus,  but  since  a  real  photon  has  rest  mass  zero,  it  can 
transfer  only  a  definite  amount  of  momentum  to  the  struck 
nucleus.   This  is  in  contrast  to  electron  scattering  where 
the  momentum  transferred  to  the  nucleus  is  limited  only  by 
the  energy  lost  to  the  nucleus. 

The  probability  that  an  incident  electron  of  total 
energy  E  is  scattered  elastically  at  an  angle  0-  (i.e.,  a, 
the  cross  section  for  elastic  scattering  at  angle  0-)  can 
be  calculated  precisely  on  the  assumption  that  the  target 
nucleus  consists  of  point  nucleons.   That  means  the  charge 
and  magnetization  densities  are  delta  functions.   The 
actual  cross  sections  differ  from  the  predicted  cross 
section  if  the  charge  and  magnetic  moment  are  spread  out, 
and  the  difference  is  usually  in  the  direction  of  smaller 
values  of  the  cross  section  since  the  smeared  out  structure 
reduces  the  field  against  which  the  particles  scatter.   The 
deviation  becomes  greater  as  the  parameter,  q,  the  momentum 
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transferred  to  the  struck  nucleon  increases.   The  (three) 
momentum  transfer  is  defined  by 

q  =  Pf  -  P±  (1-2) 

->      ->■ 
where  p.  and  p„  are  the  incident  and  final  values  of  the 

electron  momentum.   The  effect  of  structure  of  the  nucleons 

is  normally  expressed  by  the  use  of  form  factors,  G,  which 

are  believed  to  be  functions  of  the  structure  and  momentum 

transfer  only; 

p2/  2s_   cross  section  for  scattering  from  physical  particles 
q      cross  section  for  scattering  from  point  particles 

(1-3) 
Note  that  similar  ideas  occur  in  all  scattering  situations 
in  which  the  structure  of  the  scattering  center  has  an 
influence,  e.g.,  x-ray  scattering  form  factors  reflecting 
the  electronic  structure  of  atoms. 

C.   ELECTRON  SCATTERING  PROM  DEUTERIUM 

The  deuteron,  as  the  only  bound  two-nucleon  state,  has 
been  the  source  of  much  information  for  nuclear  physics, 
due  initially  to  the  great  interest  in  the  two-nucleon 
problem.   Originally  the  two-body  problem  was  thought  to 
be  of  great  importance  because  of  its  relative  mathematical 
simplicity,  and  because  of  the  hope  that  the  inter-nucleon 
force  was  additive,  so  that  knowledge  of  this  system's 
solution  would  allow  extension  to  more  complex  nuclei.   It 
is  now  generally  assumed  that  the  nuclear  force  is  not 
additive  (e.g.,  the  nuclear  force,  unlike  the  electromagnetic 
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force,  "saturates")  and  hence  the  two-nucleon  interaction  is 
not  necessarily  the  best  source  through  which  understanding 
of  the  nature  of  the  nuclear  force  can  be  gained. 

Since  free  neutrons  are  unstable  and  decay  with  a  half- 
life  of  17  minutes  to  a  proton  and  electron,  deuterium  is 
a  convenient  source  for  neutrons  as  targets.   For  this 
experiment  deuterium  in  the  form  of  solid  deuterized  poly- 
ethylene (CDp)  ,  was  used  as  a  target.   Assuming  that 
deuteron  wave  functions  are  known  from  theory,  neutron 
structure  is  inferred  through  a  comparison  of  electron- 
proton  and  electron-deuteron  scattering  results.   Conversely, 
knowing  the  nucleon  form  factors  from  theory  makes  compari- 
sons of  different  deuteron  wave  functions  possible. 

There  are  several  possible  outcomes  to  an  electron- 
deuteron  scattering  experiment: 

(1)  e  +  D  ->  e ?  +  D'      (elastic  scattering) 

(2)  e  +  D  -*  e '  +P  +  N   (inelastic  scattering,  break-up) 

(3)  Radiative  scattering 

(4)  Reactions  with  pion  production. 

Now  (2)  occurs  only  above  2.2  MeV  excitation  energy  and  is 
not  of  interest  here;  (3)  can  be  calculated;  (4)  is  not  a 
problem  for  incident  energies  smaller  than  150  MeV.   The 
theory  for  process  (1)  was  developed  initially  by  Jankus 
[3]  (1956)  in  the  first  Born  approximation,  which  should 
be  adequate  because  of  the  small  deuteron  charge.   The 
deuteron  is  described  by  nonrelativistic  wave  functions 
while  the  electron  is  treated  relativistically .   Jankus 's 
result  will  be  presented  in  Section  II,  part  B3. 
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Attempts  at  a  relativistic  deuteron  theory  were  initi- 
ated by  Jones  and  carried  forward  by  several  authors,  prin- 
cipally Gourdin,  Adler  and  Gross  [4,  5,  6,  7,  8,  9,  10,  11]. 
Results  of  their  attempts  will  be  presented  in  Section  II, 
part  B5. 

D.   THE  NEUTRON-ELECTRON  INTERACTION 

The  neutron-electron  (n-e)  interaction  is  the  term  used 
to  describe  part  of  the  electromagnetic  interaction  between 
the  neutron  and  the  electron.   This  interaction  may  not  only 
be  a  specific  interaction  between  neutron  and  electron,  but 
between  any  charged  particle  and  the  neutron.   Furthermore, 
there  exists  the  possibility  of  some  n-e  interaction  of  a 
nonelectromagnetic  (weak)  nature.   For  the  present  case 
however,  results  that  stem  from  the  interaction  of  an  elec- 
tron with  the  internal  electromagnetic  structure  of  the 
neutron  will  be  termed  "the  neutron-electron  interaction." 

Since  the  neutron  has  zero  total  charge  there  is  no 
Coulomb  interaction  between  a  neutron  and  an  electron. 
Both  particles  have  magnetic  moments  and  hence  there  exists 
a  spin-  dependent  magnetic  dipole  -  dipole  interaction 
between  them  and  also  a  velocity-dependent  interaction 
between  the  magnetic  moment  of  the  neutron  and  the  magnetic 
field  associated  with  the  convection  current  of  an  electron 
in  motion.   These  interactions  have  been  well  documented 
and  are  not  of  interest  here  [12] .   If  there  are  regions  of 
non-zero  charge-density  in  the  neutron,  (i.e.,  a  charge 
separation  exists)  then  any  charged  particle  probing  the 
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neutron  will  experience  spin  and  velocity-  independent 
electrostatic  forces.   These  are  the  forces  of  interest  in 
the  n-e  interaction. 

At  least  some  charge  separation  in  the  neutron  is 
expected  on  the  basis  that  the  free  neutron  possesses  an 
anomalous  magnetic  moment.   This  contribution  to  the  inter- 
action is  called  the  magnetic  or  Foldy  term  [13] .   The 
remainder  of  the  interaction  (due  to  charge  separation)  is 
the  result  of  the  fact  that  a  neutron  can  be  part  of  the 
time  dissociated  into  a  negative  pion  and  a  proton, 

n  «-•*■  p   +  tt"~. 

If  the  neutron  spends  80%    of  the  time  as  a  neutron  and  20% 
of  the  time  as  a  proton  +  negative  pion,  then  an  electron 
penetrating  the  'neutron'  would  see  an  electric  field 
strength  equal  to  that  produced  by  a  charge  of  0 . 2e ,  e  the 
proton  charge.   The  force  resulting  from  the  field  would 
have  a  very  short  range  since  the  tt~  would  screen  the  pos- 
itive charge  at  distances  greater  than  the  tt~  p   separation 
(The  separation  can  be  on  the  order  of  the  pion's  Compton 
wave  length  h/m  C  ~  10     cm. )   The  force  should  be 
attractive  and  can  be  represented  as  an  extremely  narrow 
potential  well.   The  purpose  of  all  the  experiments 

described  below  was  to  determine  the  well  depth  V  ,  essen- 

o 

tially  by  measuring  the  scattering  amplitude. 

The  first  of  three  techniques  to  measure  the  n-e  inter- 
action was  introduced  by  Havens,  Rabi,  and  Rainwater  [14] 
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in  19^7.   It  consisted  of  scattering  low  energy  neutrons 
(~0,1  eV)  from  liquid  lead  and  bismuth,  noting  that  the 
total  scattering  cross  section  of  an  atom  includes  an 
observable  term  arising  from  interference  between  the 
nucleus  and  the  electrons  scattering  coherently,  while  at 
higher  neutron  energies  (~10  eV)  this  interference  term 
is  almost  absent.   Results  on  liquid  bismuth  yield 

dGE  /d(q2)  p    =  0.0225  ±  0.0070  F2. 
n      q  ->0 

The  work  of  Hughes,  Harvey,  Goldberg,  and  Stafne  [151 
depends  on  the  fact  that  the  difference  in  the  refractive 
indices  for  neutrons  at  an  interface  between  bismuth  and 
liquid  oxygen  (liquid  mirror)  comes  largely  from  the  neutron- 
electron  scattering  amplitude.   The  measurement  of  the  angle 
of  total  external  reflection  gives  the  difference  between 
the  two  scattering  amplitudes  and  hence  the  n-e  interaction. 
Here  the  results  are 

dGE  /d(q2)  p    =  0.0200  ±  0.0019  F2 . 
n      q  -*-0 

The  most  recent  and  most  accurate  measurements  have 
been  by  Krohn  and  Ringo  [16]  (1966)  utilizing  a  technique 
of  Fermi  and  Marshall  [171  (19^7).   It  depends  on  the  fact 
that  in  scattering  neutrons,  the  interference  between  the 
scattering  from  the  nucleus  and  from  the  electrons  in  an 
atom  leads  to  an  asymmetric  angular  distribution  of 
scattered  neutrons  in  the  center-of-mass  system.   This 
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occurs  because  the  atomic  form  factor  depends  on  the  scat- 
tering angle.   A  summary  of  the  Krohn  and  Ringo  results  are 
shown  in  Table  I. 


TABLE  I 

Gas 

V  (eV) 
o 

dGE  /d(q2)  F2 
n 

Statistical 
weight 

Argon 

-3780  ± 250 

0.0196  ±  0.0013 

0.13 

Krypton 

-3800  ±  130 

0.0197  ±  0.0007 

0.30 

Xenon 

-3670  ±  90 

0.0190  ±  0.0005 

0.57 

Weighted 

-3720  ±  90 

0.0193  ±  0.0004 

mean 

It  is  important  to  note  that  all  measurements  agree  that 

2 
dGp  /d(q  )  is  not  zero  or  negative.   See  Appendix  A  for 

n                             dGE 
the  details  of  the  calculation  of  n   from  V  . 

d(q2) 
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II.   THEORY 

A.   THE  NONRELATIVISTIC  DEUTERON 

1.   The  Ground  State  of  the  Deuteron 

There  are  three  possible  states  of  the  two-nucleon 
system,  the  di-proton  (pp),  the  di-neutron  (nn) ,  and  the 
deuteron  (np),  of  which  only  the  deuteron  is  known  to  be 
stable.   Experimentally  measured  static  properties  are: 

(1)  Binding  Energy   W  =  2.226  ±  0.003  MeV 

(2)  Magnetic  Moment  ud  =  0.85735  ±  0.00003  nm 

where  lnm  =  1  nuclear  magneton  =  eh/2m  c  and 

p 

m  is  the  proton  mass 

(3)  Electric  Quadrupole  Moment,  Q   =  2.82  x  10~27cm2 
=  28.2  P2,  where  IF  =  10-13  cm. 

(4)  Spin  Angular  Momentum,  S  =  1ft 

(5)  Radius  rd  =  4.3  P  (see  page  23  for  definition  of  r*d) 

(6)  Mass  Md  =  2.01410  umu  (unified  mass  units) 
=  9.50  P-1 

Because  the  neutron  has  zero  charge,  the  nuclear  force 
binding  the  deuteron  cannot  be  electrical  and  neither  can 
it  be  gravitational,  for  this  force  is  far  too  weak  to 
provide  a  2.2  MeV  binding  energy.   Thus  a  nuclear  force 
that  is  strong  and  attractive  is  postulated.   By  ignoring 
some  experimentally  determined  parameters  a  simple  quantum 
mechanical  treatment  of  the  deuteron  is  possible. 

Assume  that  the  nuclear  force  is  central,  i.e.  the 
interaction  potential  between  the  neutron  and  proton  is 
some  function  V(r),  where  r  is  the  separation  of  the  two 
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nucleons.   This  assumption  is  in  conflict  with  (3)  above 
because  a  non-zero  quadrupole  moment  implies  a  nonsymmetric 
charge  distribution.   Nevertheless,  assume  the  central 
potential,  and  take  the  ground  state  of  the  deuteron  to  be 

S, .   Hence  the  Schroedinger  equation  in  the  center-of-mass 
system  is, 

2 

^  V2iK?)  +  V(r)  ip(r)  =  EiKr)  (2-1) 

where  the  reduced  mass  is, 

mm     .. 

m  =    P  "   -  |  (2-2) 

m  +  in     2 
P    n 

and  m  ,  m  are  the  masses  of  the  proton  and  neutron.   The 
p  *   n  y 

average  of  these  masses  is  M,  and  r  is  the  relative  vector 
coordinate  between  nucleons.   E  is  the  energy  of  the  rela- 
tive motion  of  the  two  particles.   Since  the  ground  state 
is  spherically  symmetric  (S  state),  i/>(r)  depends  only  on 
r  =  | r J ,  so  that 

♦  (?) -ip(r)  =  2ii2l  (2-3) 

The  quantity  u(r)  is  called  the  reduced  radial  wave  function 
Now  using  (2-3)  and  (2-2),  (2-1)  can  be  rewritten  as 

,2 
=g-  u"(r)  +  V(r)u(r)  =  Eu(r)  (2-4) 

2 
where  u"(r)  =  d  ^r)  and  E  =  -W  =  -2.226  MeV.   For  the 

dr 
square-well  potential 
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V(r)  = 


-V 


r  <  r 
—  o 


r  >  r 


(2-5) 


(2-4)  becomes 


Case  I   r  <  r 


o 


u"(r)  +  a  u(r)  =  0    (2-6a) 


where 


a  = 


■*L(V  +E)=   4(v  -W) 

h2   °     Jh2   °     (2-6b) 


The  general  solution  of  (2-6a) 
is 

u(r) =  A  sin  ar +  B  cos  ar 

(2-7) 

Requiring  that  i|>(r)  be  finite 
at  r  ■*  0  and  i/>(r)  =  0  at  r->°° 
forces  B  =  0  in  (2-7),  hence, 


u(r)  =  A  sin  ar 


(2-8) 


Case  II   r  >  r. 


u"(r)  -  y  u(r)=  0   (2-9a) 


where 


Y  = 


MW 


(2-9b) 


h< 


The  general  solution  of 
(2-9a)  is 

u(r)  =  Ce"Yr  +  Deyr    (2-10) 


the  boundary  condition  at  °° 
requires  D  =  0  so  that 


u(r)  =  Ce~yr      (2-11) 


Now  at  r  =  r  both  i>   and  ip'  (the  slope)  must  be  continuous, 
thus  (2-8)  and  (2-11)  state 

-yr 


A  sin  ar   =  Ce 
o 


o 


(2-12) 


and 


-yr 


Aacos  ar  =  -yCe 
o 


(2-13) 


The  ratio  of  (2-13)  for  (2-12)  yields 


a  cot  ar  =  y  = 

o 


MW 
u2 


(2-14) 
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This  is  a  transcendental  equation  relating  the  range  of  the 

potential  r  to  its  depth  V  .   Typical  solutions  of  (2-1*0 

for  r   vs  depth  are  given  in  Table  II. 
o 


TABLE  II 


Ranee  r  (F)  Depth  of 

Range  tq^)  Potential  (MeV) 


1.0  120 

1.5  59 

2.0  36 

2.5  25 


5.83 


The  results  of  this  rough  calculation  are  quite  good. 

Independent  measurements  of  r  yield  numbers  on  the  order 

of  1  Fermi.   If  r    2.0  F  is  accepted,  the  size  of  the 

o 

deuteron  can  be  specified  by  1/y  =  4 . 3F ,  which  is  about 
twice  that  of  the  range  r   of  the  potential.   This  explains 
why  the  deuteron  is  a  loosely  bound  system  (its  binding 
energy  is  1.113  MeV/nucleon,  vs.  8.0  MeV/nucleon  in  the 
average  nucleus). 

There  is  no  experimental  evidence  for  any  excited  state 
of  the  deuteron.   It  can  be  shown  by  simple  theoretical 
arguments  that  in  order  to  produce  a  bound  excited  state 
the  potential  well  depth  must  be  much  greater  than  in  the 
ground  state.   It  must  be  concluded  that  no  bound  state 
exists  for  £  j-    0  in  the  deuteron. 
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2 .   Review  of  Low  Energy  Neutron-Proton  Scattering 

Certain  parameters  pertaining  to  the  deuteron  have 
been  established  by  investigation  of  the  unbound  n-p  state, 
that  is  through  neutron-proton  scattering  experiments.   It 
is  clear  that  these  parameters  should  be  predicted  by  any 
good  theory  of  the  deuteron.   With  this  in  mind  this  brief 
review  is  given. 

The  n-p  scattering  cross  section  has  been  examined 
closely,  both  at  high  and  low  energies.   The  cross  section 
depends  strongly  on  the  energy  of  the  incident  neutrons. 
At  low  energies,  below  10  MeV,  the  scattering  is  essentially 
due  to  neutrons  with  angular  momentum  I   =  0  ( S-wave  scatter- 
ing) .   Thus  the  angular  distribution  of  the  scattered 
neutrons  is  isotropic  in  the  center-of-mass  system.   Since 
the  proton  has  molecular  binding  energy  of  about  0,1  eV, 
only  neutrons  with  incident  energy  greater  than  1.0  eV  are 
considered.   Hence  the  neutrons  strike  essentially  free 
protons . 

The  Schroedinger  equation  for  the  scattering  problem  is 


V2^(r)  +  ^  [e  -V(r)l  iK?)  =  0. 
fi2  L       J 


For  a  central  potential  and  for  neutrons  incident  along  the 
z-axis  the  above  becomes 


V2^(r,e-)  +  ^ 
fi2 


E  -  V(r)]  iKrGO  = 


0 


(2-15) 


which  has  the  asymptotic  solution 


2h 


<Kr,0) 


—  A 


eikz  +  f(e-)  e 


ikr 


(2-16) 


The  f(O-)  is  called  the  scattering  amplitude  representing 
the  scattered  neutrons  (waves);  Ae  l    represents  the  incom- 
ing plane  wave.   It  is  not  difficult  to  show  that 


da 
dTT 


f(e-) 


(2-17) 


Using  separation  of  variables,  ip(r,G-)  =  R(r)  Y(G-),  (2-15) 
divides  into  a  radial  and  angular  equation.   The  angular 
part  is 


-t^-k  4k   (sine-  ^  1  +  1(1+1) 


Y  =  0 


(2-18) 


with  solutions  Y„(0-)  =  P«(cos  0-)  ,  the  Legendre  polynomials 
of  order  I.      Since  the  potential  is  central  these  solutions 
are  independent  of  the  potential  V(r).   The  radial  equation 


is 


uJ(p)  + 


k2  _  u(r)   _  Ml+ii 


u£(r)  =  0 


(2-19) 


where 


,2_  2mE 
2 


,    TT,  v    2mV(r) 
and   U(r)  =  ~ — 


and  R(r)  =  u(r)/r.  Let  the  potential  V(r)  be  "turned  off." 
Then  the  solutions  of  (2-19)  are  again  the  spherical  Bessel 
functions  j-(kr).   It  can  be  shown  that  as  r  -*■  °° 


j£(kr) 


.  sin(kr-  £  it/2)   „v   f „ 
kr 
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and  hence 

UZM    r -v  a,  "  sin(kr  -  £tt/2).  (2-20) 

When  V  is  "turned  on"  but  r  is  still  greater  than  r  ,  the 
solution  of  (2-19)  must  still  correspond  to  that  of  a  free 
particle,  and  hence  be  the  same  as  (2-20)  except  for 
possibly  a  phase  factor.   Thus 

Vr)  r  ->  co-  si"(kr  -  4r  +  <V    "V  on".       (2-21) 
The  general  solution  to  (2-19)  for  "V  on"  at  r  >  r   is 

CO 

♦  (r,*)    =    £     C,   sln(kr-£ir/2+fi£)    p£(cos0) 

£  =  0 

(2-22) 

CO 

2,   B0Ppi;cose-) . 

£=0 

But  this  must  be  the  same  as  (2-16).  Now  using  a  Rayleigh 
expansion  on  e  in  (2-16),  and  writing  all  quantities  as 
complex  exponentials,  equating  (2-16)  and  (2-22)  yields 

r  1*(2&+I)e   : 


and 

f(G-)  =  i^  (2£  +  l)e  l    sin6JlPjl(cose-).  (2-23) 

I 
The  square  of  the  scattering  amplitude  is  the  differential 

scattering  cross  section  (2-17),  hence 

cT  =  f  m^1   sinO-  dO-  =  ?|^  (2£  +  l)sin26£.      (2-24) 
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Note  that  S-wave  scattering  (£  =  0)  is  predominant  for 
energies  less  than  10  MeV.   By  taking  the  range  of  nucleon 
force  to  be  2F  and  assuming  a  0.5  MeV  incident  neutron 
(in  center-of-mass  frame  or  1  MeV  in  the  lab)  it  can  be 
seen  by  forming  the  ratio  of  the  coefficients  B.  in  (2-22) 
that 

2 


Bi 


B 
o 


=  (kr  )2 


This  states  that  only  about  9%   of  the  scattering  is  due  to 
neutrons  with  £  =  1. 

In  (2-24)  let  I   =    0,  then 

)i^     n       /sin  6  \2 

As  k  -*■  0.  sin  6   approaches  6  .   In  order  that  am  remain 

o  o  T 

finite  6  /k  must  approach  some  constant,  called  -a  by 
convention.   Thus 

aT  =  ^Tra2  (2-26) 

where  a  is  called  the  scattering  length.   Consideration  of 
the  solution 

6 

u  (r)  ~  sin(kr  +  6    )    =  sink(r+  -r— )  .  ^  ■  k(r-a) 
o  o  k   k-*-  0 

demonstrates  the  fact  that  the  asymptotic  solution  becomes 
a  straight  line,  and  that  the  scattering  length,  a,  is  the 
intercept  on  the  r-axis  and  is  obtained  by  extrapolating 
the  radial  wave  function  inside  the  potential,  well  beyond 
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the  range  of  force  r  .   Figure  1  illustrates  the  signifi- 
cance of  the  sign  of  the  scattering  length.   The  scattering 
length  a  is  positive  if  the  state  is  a  bound  state. 


(a) 


Figure  1 


(b) 


Positive  scattering  length 
indicates  a  bound  state 


Negative  scattering  length 
indicates  an  unbound  state 


The  energy  dependence  of  the  scattering  cross  section 
can  be  expressed  in  terms  of  the  scattering  length  a,  and 
another  parameter  which  has  the  dimensions  of  length  and 
is  known  as  the  effective  range.   Denote  this  parameter 
by  p,  where  p  is  of  the  same  magnitude  as  the  range  of 
nuclear  force  r  .   Effective  range  theory  predicts  the 
phase  shifts  6.  as  a  function  of  energy.   Consider  S-wave 
scattering  only.   The  wave  equations  for  two  energies 
E  =  h2k2/2m  and  E  =  0  are 
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Case    I      V  i-    0 


u"+  (k2-^)u  =  0  (2-27a) 


u"  -  M-  =  0 


(2-27b) 


fi' 


multiply    (27a)    by    u0 ,    (27b)    by 
u  and   subtract,    getting 


h   (uuo-uoU,)  °  k'uuo   (2~28) 


Case  II   V  =  0 


v"  +  k  v  =  0       (2-29a) 


v^  =  0  (2-29b) 


multiply  (29a)  by  vQ ,  (29b) 
by  v  and  subtract,  getting 


-r-(vv'  -  v  v')  =  k2v  v  (2-30) 
dr   o   o        o 


Subtract    and   integrate    the    above   equations.      Then 

p=oo 


(  uu '  -u   u '  -  vv '   +v   V  '  ) 
o        o  o        o 


r=0 


/ 


=   k      /    (uu   -vv    )dr.    (2-31) 
o        o 

"0 


Because  u  =  v,  u   =  v   for  r  >  r  ,  the  expression  on  the 

o    o  o ' 

left  hand  side  contributes  nothing  at  the  upper  limit. 

Furthermore  u(0)  =  u  (0)  =  0,  hence  (2-31)  is 


(vv'  -  v  v') 
o    o 


=  k1 


r=0 


oo 

/(uu  -vv 
o   o 
0 


)dr. 


(2-32) 


In  analogy  with  (2-20)  the  solutions  u  and  v  can  be  written 

sin(kr+6  ) 


u(r) ^v(r)  =  C  sin(kr+6  )  =  — ^  .  R 

1    -p  •*■  °°  o       sm  6 


where  C  is  a  normalization  constant.   Now  (2-32)  becomes 


-  —  -  k  cot  6   =  k  I     (uu  -vv  )dr 


i 


and  in  the  limit  of  zero  energy  neutrons  u  -  u  ,  v  -  vq  so 


that 


k  cot  6  ="—-   +  p  — 
a      2 
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where 


00 

/ 


p  =  2     (v2  -  u2)dr 
o    o 


0 

is  the  effective  range. 
Hence  the  cross  section  is  given  by 

°t  =  ~ — ~rn  (2-33) 

k2  +  (i-£|-) 

2        2 
Besides  being  a  function  of  k   =  2mE/h  ,  the  cross  section 

is  expressed  in  terms  of  the  two  parameters,  effective  range 

p  and  scattering  length  a.   The  cross  section  is  independent 

of  the . form  and  shape  of  the  potential.   Both  p  and  a  can 

be  determined  by  measuring  the  cross  section  at  different 

low  energies  of  the  incident  neutron. 

The  theoretical  estimate  of  the  value  of  the  scattering 

cross  section  for  I   =  0  is  far  below  the  experimentally 

observed  value.   In  fact  for  a  square  well  of  reasonable 

width  and  depth 

oT  -  4.5  x  10~2Z|  cm2 

while 

a      .      ,  -  20.4  x  10    cm  . 
experimental 

This  conflict  was  removed  when  Wigner  realized  that  the 
theory  for  aT  is  based  on  a  potential  whose  "strength"  is 
calculated  from  the  deuteron,  where  the  spin  is  1,  i.e. 
neutron  and  proton  are  parallel }  while  scattering  experi- 
ments are  done  using  randomly  oriented  nuclei.   The 
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scattering  experiments  then  represent  an  average  over 
parallel  (triplet)  and  anti-parallel  (singlet)  states.   The 
weight  factor  for  the  singlet  and  triplet  cases  are  1/k    and 
3/*J  respectively,  hence  the  total  cross  section  for  scatter- 
ing is 


%  =  ¥  °t  +  ¥ 


a 


s 


where  a,  and  a   are  the  triplet  and  singlet  cross  sections. 

Using  (2-33)  and  working  backwards  assuming  a,  =  20.^  barns, 

1/Y+.  =  ^.3  P,  and  r,  ^  r   -  2  F  it  can  be  determined  that 

W  =  0.066  MeV  with  1/y   =  25  F,  where  the  scattering  length 
s  s 

a   is  negative.   Hence  the  singlet  state  is  not  a  bound 
deuteron  state. 

3 .   Pheonomenological  Two-Nucleon  Potentials 

On  the  evidence  presented  so  far  it  can  be  concluded 
that:   (1)  The  deuteron  is  predominantly  in  a  S-.    state 
since  the  solutions  to  the  Schroedinger  equation  for  a 
central  potential  are  quite  good;  (2)  since  \i-.^)i      +  y  , 
and  Q ,  j-    0 ,  there  is  reason  to  guess  that  noncentral  forces 
are  acting  (the  inequality  of  magnetic  moments  in  a  complex 
problem  and  could  well  have  other  interpretations);  (3)  the 
forces  between  neutron  and  proton  are  different  in  the 
singlet  and  triplet  states,  implying  that  nuclear  forces 
are  spin-dependent.   There  is  a  point  of  confusion  on  (3). 
Spin-dependent  does  not  imply  that  this  contribution  is  not 
central  in  nature.   It  means  that  the  usual  central  scalar 
potential  is  different  according  to  the  relative  orientation 
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of  spins,  a  situation  that  has  no  direct  analogy  either  in 
electromagnetism  or  atomic  physics. 

A  realistic  deuteron  potential  must  incorporate  these 
ideas  and  more.   The  potential  must  be  symmetric  in  the 
coordinates  of  the  two  nucleons  and  invariant  under  rota- 
tions and  reflection  of  their  space  coordinates.   This 
means  the  potential  is  a  scalar. 

Some  of  the  more  detailed  features  of  the  nuclear 
potential  can  be  directly  related  to  the  spin  and  isospin 
dependence  of  the  nucleon  interaction.   An  important 
contribution  for  the  nuclear  symmetry  potential  comes  from 
the  exchange  character  of  the  nuclear  forces .   Thus  a 
Serber  exchange  potential,  Pr,  for  the  central  interaction 
would  be 

Vc  =  |  (1  +  Pr)  V(r)   or 

Vc  =  |  [1  -  |  (1+  t1-t2)(1  +a1'a2)]V(r)        (2-34) 

where  t  ,  t_  and  a, ,  a?  are  the  isospin  and  Pauli  spin 
matrices  respectively. 

The  occurrence  of  a  rather  strong  spin-orbit  force  in 
the  nuclear  interaction  gives  rise  to  a  spin-orbit  coupling 
in  the  two  body  force; 

vls  =  irvLs(r)(W  xcPi-^-cVV-     (2-35) 

This  VT Q  contribution  is  based  on  Wigner's  idea  that 
velocity-dependent  forces  are  acceptable  when  they  depend 
only  on  the  relative  momentum  p  of  the  two  nucleons. 
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All  these  facts  taken  together  yield  a  general  potential 
V  =  Vc(r)  +  Va(r)  o1-a2   +  VT(r )S12  + VLS (r)L- S  (2-36) 

where  a-\'c?   represents  the  spin-spin  interaction  and  S-«  is 
the  tensor  operator 

S12  =  3(a1-r)(a2-r)  -  (a^o^.  (2-37) 

The  form  of  the  tensor  operator  can  be  inferred  by  analogy 
with  the  classical  problem  of  the  potential  energy  of  an 
electric  dipole  d,  in  the  field  of  another  dipole  dp  at  a 
distance  r.    Here 


r  3  h(d1-r)(d2-?)  -  d-^dp].  (2-38) 


The  r    comes  from  the  fact  that  the  electric  charges 
interact  with  a  potential  proportional  to  r~  .   But  the 
r-dependence  of  "mesic  charges"  which  are  thought  to  be  the 
source  of  nuclear  fields  is  unknown.   Hence  the  r   depen- 
dence must  be  replaced  by  an  arbitrary  function  of  r,  VVp(r) 

The  structure  of  the  nuclear  force  is  closely  related 
to  the  properties  and  interactions  of  the  whole  family  of 
strongly  interacting  particles.  It  has  not  been  possible 
so  far  to  derive  the  forces  between  these  particles,  or 
their  masses,  from  simple  assumptions  regarding  the  basic 
structure  of  the  strong  interaction.  Nevertheless,  certain 
relations  between  the  interactions  and  masses  of  the 


s12. 


*See  Appendix  B  for  properties  of  the  tensor  operator 


33 


strongly  interacting  particles  can  be  established.   Of 
particular  significance  is  the  relationship  between  the 
nuclear  force  at  large  distances  and  the  pion-nucleon 
interaction. 

The  force  at  large  distances  can  be  described  in  terms 
of  an  exchange  of  mesons,  in  a  similar  manner  as  the  elec- 
tromagnetic interaction  can  be  analyzed  in  terms  of  photon 
exchange.   The  special  role  of  pion  exchange  is  due  to  the 
smallness  of  the  pion  mass  as  compared  with  that  of  other 
strongly  interacting  particles.   The  interaction  associated 
with  the  exchange  of  a  particle  of  mass  m   is  limited  to 
a  range  of  the  order  of  the  compton  wave  length  X  =fi/m  c. 

In  fact,  the  intermediate  states  involved  in  such  an 

2 
exchange  have  an  energy  of  at  least  m  c   and  so  are  limited 

2 
to  a  duration  of  order  fi/m  c  .   During  this  time,  the 

emitted  particle  cannot  travel  farther  than  h/m  c  and  thus, 

the  interaction  is  expected  to  decrease  strongly  for 

distances  greater  than  X  .   If  n  particles  are  exchanged 

simultaneously,  the  corresponding  intermediate  energies 

2 

are  n  x  m  c   and  the  range  is  X  /n.   The  interaction  at  the 
7T  °        c 

largest  distances  is  thus  determined  by  the  exchange  of 
single  mesons. 

The  asymptotic  form  of  the  one-pion  exchange  potential 
(OPEP)  is  uniquely  given  by  the  mass  and  symmetry  properties 
of  the  pion  (J*  =  0~,  T  =  1) 


3^ 


VOPEP   =    3   h^"  m7TC    (T1    T2} 


(o1-o2)   +  S12(l  +  |  +  -|-) 


-X 

~-   (2-39) 


where  x  =  m^  c/ft  =  0.70  P  -1,  f2/hc  =  0.081  ±  0.002. 

The  strength  of  the  potential  is  determined  by  the  coupling 
constant  f  for  the  process  N  ->  N  +  it  .   The  meson  field  is 

IT      — 

pseudoscalar  since  J   =  0  ,  and  isovector  since  T  =  1. 

Equation  (2-39)  reduces  to  (2-36)  for  the  deuteron. 

Early  definitive  work  on  the  phenomenological  nucleon- 

nucleon  potential  in  the  OPEP  approximation  was  carried 

out  by  Glendenning  and  Kramer  [18].   They  attempted  to 

construct  triplet  (spin)  -  even  (parity)  potentials  that 

are  asymptotic  to  the  OPEP  and  are  modified  in  the  inner 

regions  so  that  the  deuteron  properties  are  obtained.   The 

results  of  an  8-parameter  fit  to  an  equation  similar  to 

(2r-39)  are  found  in  Table  III,  entries  1-9.   Coefficients 

in  their  potential  are  varied  such  that  the  empirical  values 

of  Qd  =  28.2  P2  and  a,  =  5 .kk    ±  0.02  F  are  reproduced.   A 

hard  core  radius  r  was  chosen  so  that  the  potential  has  a 

c  ^ 

bound  state  at  W  =  2.226  MeV.  (The  hard  core  is  discussed  on 
page  3  8). 

Many  other  authors  have  contributed  to  the  data  [19,  20, 
21] .   One  of  the  most  useful  has  been  the  Hamada-Johnston  (HJ) 
[22]  potential  given  by 

V  =  Vc(r)  +  VT(r)S12  +  V^L-S  +  VLL(r)L12      (2-40) 
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TABLE  III 


Potential 
Identification 

re(f) 

at(F_1) 

P(F) 

Qd(P2) 

pDm 

X 

GK1 

0.4815 

5.376 

1.717 

28.79 

7.415 

0.02714 

2 

0.2466 

5-368 

1.709 

28.21 

6.281 

0.02693 

3 

0.47^7 

5.456 

1.812 

28.17 

5.970 

0.02595 

4 

0.4425 

5.364 

1.703 

28.22 

7.103 

0.02677 

5 

0.3632 

5.366 

1.705 

28.02 

6.710 

0.02663 

6 

0.3924 

5.384 

1.726 

28.12 

5.957 

0.02676 

7 

0.5369 

5.477 

1.836 

28.14 

6.029 

0.02564 

8* 

0.5007 

5.413 

1.760 

28.18 

5.622 

0.02654 

9* 

0.4329 

5.373 

1.715 

28.29 

7.425 

0.02670 

HJ 

0.485 

- 

1.77 

28.50 

6.97 

0.02656 

FL 

0.735 

- 

- 

- 

4.60 

- 

Experiment 

- 

5.40  + 

- 

+  + 
28.20'T 

- 

- 

X  is  the  asymptotic  D  state/S  state  ratio 

*  GK8  and  9  include  the  spin-orbit  term  in  (2-34)  while  GK1-7 
are  without  it . 

t   N.  K.  Glendenning  and  G.  Kramer,  Phys.  Rev.  126,  2159,  1962, 

tt  G.  P.  Auffray,  Phys.  Rev.  Letters  6,  120,  196I. 
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where 

L12  =  (o^'c^)  L^  -i[(a1-L)(a2-L)  +  (a2-L)(a1«L) 

=  (5LJ+  0^0^    l/  -  (L-ST 

and  contains  central  V  .  tensor  Vm,  spin-orbit  VT c .  and 
second-order  spin-orbit  VT T  components.   It  is  from  this 
potential  that  the  later  used  Partovi  wave functions  have 
been  derived.   The  radial  functions  are  restricted  by  the 
condition  that  at  large  distance  the  central  and  tensor 
potentials  should  be  described  by  the  OPEP .   The  radial 
forms  are : 

Vc  =  v0(T1-T2)(a1-ff2)  Y(x)[l  +  acY(x)  +bcY2(x)] 
VT  =  vo(T1'T2)(a1'a2)  Z(x) [1  +aTY(x)  +bTY2(x)] 

VLS   =    ^LS    V2(x)[1+bLS   Y(X)] 

VLL   "    SLL   Vo   ^    [l+aLLY(x)  +  bLLY2(x)l 
x 

v      =   i  J-  m   c2    =    3.65   MeV 
o        3   he      it 

m   c  -x 

tt  r  „/    x         e 

x  =  it  r  =  itw  •   Y(x)  =  — 

z(x)  =  1  +  3.  +  _L  . 

x        x2 

In  addition  the  potential  has  been  assumed  to  have  a  com- 
ponent giving  rise  to  an  infinite  repulsion  at  radius 

x   =  0.3^3  (r   =  0.485  F) . 
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The  values  of  the  quantities  V  ,  V™ ,  etc.  have  been 
calculated  for  four  different  types  of  potentials.  The 
results  are  shown  in  Table  IV. 


TABLE  IV 

Vc(MeV)    VT(MeV)  VLS(MeV)  VLL(MeV) 

Singlet  even     -1460      -  -         -k2 

Triplet  even     -  207     -6^2        3^        668 

Singlet  odd      2371      -  -  -6688 

Triplet  odd     -   23      173  -1570  -1087 


Values  of  HJ  potential  at  r  =  0 . M85  F. 

c 


The  "hard  core"  mentioned  above  is  a  phenomena  that 
cannot  be  investigated  through  low  energy  nucleon-nucleon 
scattering.   In  the  neighborhood  near  and  above  200  MeV 
repulsive  components  in  the  nucleon-nucleon  interaction  are 
seen.   One  way  of  introducing  a  repulsive  interaction  (which 
might  account  for  nucleon  saturation)  is  in  terms  of  a  very 
strong  repulsive  potential  of  very  short  range  -  "the  hard 
core."   The  infinite  short  ranged  repulsion  is  only  one, 
rather  extreme,  way  of  accounting  for  the  observed  effects. 
Finite  repulsive  potentials,  "soft  cores,"  as  well  as  inter- 
actions depending  explicitly  on  relative  velocity  can  account 
for  some  experimental  results. 
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A  fairly  recent  idea  by  Feshbach  and  Lomon  (FL)  [231  in 
1967,  is  to  replace  the  interior  region  by  appropriate 
boundary  conditions  at  some  finite  radius.   As  long  as  the 
interactions  responsible  for  these  boundary  conditions 
involve  energies  much  greater  than  the  bombarding  energy, 
it  may  be  expected  that  the  boundary  conditions  will  be 
approximately  energy  independent.   The  FL  model  also  employs 
the  (TPEP)  two  pion  exchange  potential. 

4 .   Solutions  of  the  Non-Central  Force  Deuteron  Problem 
Solutions  to  the  non-central  force  deuteron  problem 
are  obtained  by  substituting  the  potential  v"npFp  into  the 
Schroedinger  equation  and  solving  for  the  wave  functions. 
Angular  momentum  and  parity  arguments  can  determine  the 
deuteron  ground  state  angular  wave  function. 

The  total  angular  momentum  of  the  deuteron  is  J  =  1 
and  it  is  in  a  definite  parity  state.   The  value  J  =  1  can 
be  obtained  from  different  combinations  of  the  orbital 
angular  momentum  I,    and  the  spin  s  such  as: 

_3c 


£  =  0    s  =  1- 

•  s  =  0 
I   =  1  , 

.  s  =  1- 

I   =  2    s  =  1 


^P. 


~3, 


~3D. 


and  no  others.   If  the  deuteron  is  a  mixture  of  these  it 

■3-5      1     ? 
can  only  be   S,  +   D,  or   P,  +   P,  since  the  combination 

must  have  a  definite  parity.   Since  an  almost  spherically 
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symmetric  ground  state  exists,  the  obvious  choice  is  the 


S  state  combination.   It  is  not  difficult  to  show  that  the 


1     3 
odd  parity  states  (  P.,  +   Pn  )  yield  the  wrong  values  for 


Qd  and  yd« 


M 


The  angular  momentum  wave  function  14   ^-j    for  a  two 
particle  system  in  the  orbital  angular  momentum  state  i   and 
spin  state  S,  forming  the  total  angular  momentum  J,  with 
Z^-component  M  is  given  by 
M 


m0   M 
£    s 


£SJ  "  I  <*Sm,Ms|jM>  Y£   X 

Ms>m£ 
M 


(2-Ml) 


M 


*  s 

where  the  Y„   are  the  spherical  harmonics,  x    are  the  spin 

wave  functions  and  <&  Sm„  M  I JM>  are  the  Clebsch-Gordon 

I      s 

coefficients.   Two  spin  h   particles  have  the  following  spin 
states : 

Triplet  (s  =  1) 

X^  =  <%¥s%|ll>a(l)a(2)  =  a(l)a(2)  (2-42a) 


X?  =  <hh~hh\ 10>a( 1)3(2 )+ <hhh~h\l0>a(2) 6(1) 


_1_ 
/2 


a(l)3(2)  +  a(2)  6(1) 


X71=  <^35-35-35|i-l>B(i)3(2)  =  6(1)3(2) 


(2-H2b) 
(2-42c) 


Singlet  (S  =  0) 


Xq  =  <¥^-1-2|00>a(l)6(2)  +  <¥r-V^|00>3(l)a(2) 


0    J^_ 
X°    /2 


a(l)3(2)  -  3(l)a(2) 


(2-43) 
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where  a  =  (_)  and  $=(-.)  are  "spin  up"  (proton)  and  "spin 
down"  (neutron)  states  of  the  nucleon.   The  arguments  of 
the  a's  and  B's  stand  for  particle  number  one  or  two.   The 


angular  momentum  wave  functions  IL 
under  consideration  are 


M 


for   the    two   states 


^r      J  =    1,      M  =    1,      I  =    0,      s    =    l 


01J   =    <010l|ll>    Y°a(l)a(2)    =   Y°a(l)a(2)  (2-44) 


3D1:      J   =    1,      M  =    1,      £=2,      s=l 


1-M        M 
s        s 


J    =^     <21(l-Ms)Ms|ll>    Y2      S    x 


M 


=  ^JY2    PCD3C2)    --J5  Y^    -i-[o(l)B(2)+3(l)o(2) 


+  \1^  Y2    a(1)a^2)- 


(2-45) 


is 


The   Schroedinger  equation   in   the    center-of-mass    system 


where 


,2      p 

=2-  V^  +  V     +  vmsno 
2m  c  T    12 


(J)      =   E<J> 


(2-46) 


<j)    =    (j>      +   <f> 


_      u ( r ) o /    1      +   w(r) 


y1  ?y     ^ 
Oil  r       7  211 


(2-47) 


is   the   ground   state    deuteron  wave    function.'      For   ease   the 
linear   spin-orbit    and   quadratic    spin-orbit    parts    of   the 
potential    (2-40)    have   been   ignored.      Rearranging    (2-46) 
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#v2  +  V 
2m  c 


u(r )  Ojj       1        w(r) 


Oil 


y21x 


+   V, 


T 


y2ii  +  ?K7oii-^2ii 


F^ii  +  F 


=    0 


(2-48) 


In  writing    (2-48)    the    following  properties    of   S, p   have   been 
used : 


S12^011         ^ 


(2-49a) 


and 


y2ii-^?c 


si2y2n=  ^foii 


1 

211 


(2-ll9b) 


Note   that 


u2   =  _L  A   r    2   -A>  1/ 

2    3r    lr      9r;    ~  <2    2 
r        L  "fir 


(2-50) 


where 


2   _      -h' 


IT    = 


sin20- 


sin   q.  ^   (sin   Q.    *_)    +     9 

dtp 


(2-51) 


and 


p   mp  m 

i/Y^    (0-}(f))    =    *U  +  l)fi      Y£*    (e-,4>). 


(2-52) 


Hence 


v   VTon; 


1    dfu 

r   a    2 
dr 


1 
Oil 


(2-53) 


and 


^1%A> 


(l    dfw        1    2(2+1) 

I  r    ,    2         r  2 

\      dr  r 


w 


1 
211 


(2-54) 
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Multiply  (2-48)  by  U   0ii>  integrate  over  the  solid  angle, 


yM 


£SJ 
to  get 


-f>2 

-—■  u"(r)  +  (Vc  -  E)u(r)  =  -/8  VT  w(r)  .        (2-55) 

Similarly,  multiply  (2-48)  by  l{    p,,  and  integrate  to  get 


2m 


[wn(r)  _  6  wi|lj  +  |- 


V   -  2Vm  -  E 
c      T 


w(r) 


=  -  /8  VTu(r).  (2-56) 

These  coupled  differential  equations  were  first  obtained  by 
Rarita  and  Schwinger  [24] . 

A  method  for  solving  these  equations  has  been  described 
by  Hamada  [21]  ,  who  states  that  (2-55)  and  (2-56)  have  two 
sets  of  linearly  independent  solutions.   Thus  choose  (u.,  ,w.  ) 
and  (Up,Wp)  with  the  asymptotic  requirements  that  the  first 
solution  behave  as  an  S  wave  function  and  the  other  as  a  D 
wave  function: 


-ar 
e 


(2-57a) 


and     /u0\       /  0 

e"ar[l  +  3/ar  +  3/(ar)2]  j. 


(2-57b) 


Starting  with  the  behavior  at  large  r,  integrate  (2-48) 
inward  to  get  the  two  solutions  for  arbitrary  r.   Then  u(r) 
and  w(r)  are  obtained  as  a  linear  combination  of  these  two 
solutions : 
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u(r)1 


w 


Ml     d+x2)h 


U1(T)VJ   U2(r) 


w1(r) 


Wp(r) 


(2-58) 


The  requirement  that  u  and  w  in  (2-58)  vanish  at  the  core 
radius  r   gives  an  equation  for  r  which  is  solved  by  an 
interpolation  procedure.   The  asymptotic  D  state/S  state 
ratio  X  is  then  determined  from 


A  =  - 


Mr  J 


Wn  (r  ) 


u2(rc}      w2(rc}  " 


(2-59) 


The  solutions  (2-58)  are  found  in  practice  by  numerical 
methods.   Now  having  the  wave  functions  in  hand  the  magnetic 
moment  and  quadrupole  moment  can  be  predicted  by  taking  the 
expectation  values 


/■ 


<yz>  =  j  <T  yz<(>  dr, 


y  =  —  I      +gs    +gs 
z   2   z    bp  pz    ^nnz 


(2-60) 
(2-61) 


and 


■/■ 


<Qd>  -  /  <T  Qd  <f>  dr, 


A  P  P 

Qd  =  3Z^  -  r  . 


(2-62) 
(2-63) 


Typical  results  are  listed  in  Table  III  for  various  poten- 
tials.  (See  Appendix  C  for  example  calculations  of  <y  > 
and  <Qd>. ) 

B.   THE  SCATTERING  OF  ELECTRONS  FROM  PROTONS  AND  DEUTERONS 

Traditionally  electron  scattering  cross  sections  have 
been  handled  through  the  techniques  of  first  order 


Hk 


time-dependent  perturbation  theory,  the  main  result  of  which 
is  the  Fermi  Golden  Rule  No.  2.   This  rule  relates  the 
transition  probability  per  unit  time,  r,  to  the  square  of 
the  matrix  element  for  the  transition.   Specifically, 

r  =  2*  |<f|Hl|i>|2  §  (2-64) 

where  dN/dE  is  the  density  of  states  like  | f >  around  the 
energy  E,  and  H-j-  is  the  interaction  Hamiltonian.   The  cross 
section  is  then  defined  as 

Y 
incident  particle  flux"  ^  ~  ^ ' 

Quantum  electrodynamics  (QED),  which  is  the  interaction  of 
the  Dirac  spinor  field  with  the  electromagnetic  field,  pro- 
vides the  mechanism  for  computing  the  cross  section  via 
S-matrix  theory.   In  general  terms  the  S-matrix  is  a  unitary 
operator  that  transforms  the  initial  state  |i>  into  the 
final  state  | f > ,  i.e.,  S|i>  =  | f > ,  where 


sfi 


oo 

=  y  \   (-i)n  /dV  ..  ./d^x  t[ht(x,)...Ht(x  ) 
Z-rnI       /    1    /    n  I  I   1     In. 

n=0  (2-66) 


and  x  is  the  chronological  operator  that  orders  the  inter- 
action with  respect  to  causality.   The  expression  (2-66)  is 
integrated  over  all  four  space  and  hence  is  not  physically 
measurable.   If  Sf.  is  normalized  to  unit  space-time  volume 
then  it  can  be  written 


3   I 

(P,  -  P,)K,  |  (2-67) 


VT        v^f   *i''  fi 
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2j 
where  6   is  a  four  dimensional  delta  function  expressing  the 

conservation  of  four  momentum  in  the  transition  from  |i>  to 

| f > .   The  cross  section  for  a  transition  from  a  two-particle 

initial  state  |i>  into  some  group  of  final  states  is  in  the 

cent er-of -mass 


1        V     f 


'■^i^-r-f-r-  2      J    ^(pf-Pi)|Mflrdpr..dpn 

t? ("*   ^      t?  r "**    ^    final    final 

Mpl;      Mp2;    spins   momenta  (2-68) 

( Integration  over  suitable  variables  yields  the  desired  dif- 
ferential cross  section  da/dft. 

At  the  heart  of  (2-68)  is  the  matrix  element  |Mf.|.   It 
is  here  the  character  of  the  interaction  is  found  while  the 
other  factors  represent  kinematical  quantities.   The  form 
of  Mf.  can  be  surmised  by  examination  of  the  Feynman  diagrams 
for  the  process  under  consideration. 

Much  of  QED  is  conveniently  described  by  Feynman  diagrams. 
Here  the  diagrams  will  be  mainly  used  to  supplement  the 
understanding  of  the  electromagnetic  interaction,  but  their 
use  is  certainly  not  limited  to  this  role.   If  the  diagram 
can  be  drawn,  then  there  exist  specific  rules  for  construct- 
ing the  S-matrix  therefrom. 

Electron  scattering  is  graphically  represented  by 
Figure  2.   Both  graphs  are  equivalent,  but  (a)  is  most  often 
seen  since  it  is  just  a  little  more  suggestive  of  the 
scattering  process  in  the  center-of-mass  frame. 


k6 


A 


(a) 


(b) 


Figure  2 


The  time  axis  is  directed  upwards  unless  otherwise  noted. 
The  horizontal  axis  can  represent  a  space  or  momentum  inter- 
val, the  distinction  between  the  two  being  of  no  importance 
here.   The  arrows  on  the  lines  point  in  the  direction  of 
increasing  time  for  a  particle.   If  the  arrows  were  reversed 
the  particle  then  represented  would  be  the  antiparticle  of 
the  original  one.   Figure  3  sometimes  is  said  to  represent 
particles  in  the  "scattering  channel"  while  the  "crossed  or 
annihilation  channel"  is  represented  by  the  same  diagram 
rotated  by  tt/2  . 


rotated  ccw 
by 

Z 


scattering 
channel 


Figure  3 


annihilation 
channel 
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Lines  beginning  and  ending  within  a  graph  represent  "virtual 
particles."   Table  V  below  lists  some  of  the  more  important 
graph  elements  seen. 

TABLE  V 


Graph 
Element 


Factor  in  the 
S-Matrix 


Interaction 


gr-  u(p)  ,  /gr-  u(p 


P 


P 


annihilation  (creation) 
of  an  electron 


m   /  x   /  m  — 

E"  v(p)>   E~  Vp 
P      V  P 


annihilation  (creation) 
of  an  antielectron 


e  (k) 


V? 


10 


k 


annihilation  (creation) 
of  a  photon 


iyp  -  m 

2  ,   2 
p  +  m 


electron  progator 


photon  progator 


i  2  .  2 

k  +  m 


meson  progator 


Yy<$(p-p'-k) 


VnV 


YM0M6(p-p'-k) 


photon-electron  vertex 


photon-nucleon  vertex. 
Bubble  or  cross  hatched 
circle  represents  ignor- 
ance of  the  true  inter- 
action; mathematically 
this  ignorance  is  repre- 
sented by  form  factors 
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1 .   Elastic  Electron-Proton  Scattering  from  Point  Protons 
To  calculate  the  elastic  electron-proton  (e-p)  scat- 
tering cross  section  from  point  (Dirac)  protons  in  the  first 
Bern  approximation  it  is  assumed  that  (1)  the  electron  is  a 
physical  point,  and  (2)  the  one  photon  exchange  process  des- 
cribes elastic  scattering.   The  graph  of  Figure  2  is  then  valid. 


Figure  k 

From  the  figure  define  the  kinematical  quantities  p.,pf, 
P*»Pf  as  the  four-momentum  of  the  incident  and  final 
electron  and  proton  respectively;  q  is  then  the  four- 
momentum  transfer 


q2  5  (p.  -  pf)2  =  (P.  -  Pf)2. 


(2-69) 


2  * 

Choose  the  convention  that  q   <  0  for  electron  scattering. 

By  looking  at  Figure  H    it  can  be  seen  that 


The  metric  used  is  q   =  q   -  q  .   In  the  lab  frame  the 
recoil  energy  given  the  struck  particle  is  small^with  respect 
to  the  total  incident  electron  energy  thus  q^ -  -q2 <  0.   The 
four  vector  q2  is  then  said  to  be  spacelike.   In  the  annih- 
ilation channel  q2  >_  4M2_^>  0  and  is  timelike.   Another 
metric  commonly  used  is  q2=q2-q2.   The  different  metrics 
yield  cross  sections  and  form  factors  with  some  sign  dif- 
ferences.  This  possibly  ambiguity  is  being  remedied  by  use 
of  the  variable  t  for  ±q2  in  many  newer  papers. 


^9 


sfi  ~  Mfi  =  k  4  jP  (2-7o) 

q 
which  represents  the  interaction  of  the  electron  four- 
current  j   with  the  four-current  of  the  proton  J   via 
virtual  photon  progator  —  2. 

Equation  (2-70)  is  the  equivalent  of  taking  the  inter- 
action Hamiltonian  to  be 

HjCx)  =  J  (x)  AP(x)  (2-71) 

where  A  (x)  is  the  electromagnetic  field  of  the  proton, 
which  in  turn  depends  on  the  proton  current  J  (x).   Regard- 
less of  the  point  of  view  the  cross  section  (2-68)  becomes 

a  -  6lj(pf-  p.+  Pf  -Pi)|j   S"  jlJ|2'  (2-72) 

q 

QED  predicts  the  form  of  the  current  for  the  electron, 

Jy ►u(pf)Yyu(p.)  (2-73) 

and  by  analogy  for  the  proton 

Jy -ur(Pf)yyU(P1)  (2-74) 

where  u,U  and  u,U  are  the  incident  and  final  electron  and 
proton  spinors ,  respectively;  and  Y  is  a  4  x^  matrix.   Sub- 
stitution of  the  currents  into  (2-72)  and  carrying  out  the 

i    1 2 
trace  calculation  of  |Mf.|   followed  by  a  sum  over  final 

spins  and  an  integration  over  final  momenta  will  finally 

yield  the  result 
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*°\  J*!)2   2°44(l  +  -^  sin2   tWl-fifal  tan2 


(2-75) 


p 

Proton 


where 


2  \2     2  —  /   \ 

If  I  C°Sj.  I  =  Mott  cross  section  =  [§£]  (2-76) 

2E/  sin11  £  \dfi/Mott 


and  [1  +  ^p  sin2  |\   =  recoil  factor  =  £  (2-77) 

are  the  results  of  the  kinematical  quantities  in  (2-68). 
The  square  of  the  invariant  matrix  element  Im^.I  yields  for 
a  Dirac  proton 

2  2 
|Mf± |2  =  1  -  ^-^2  tan2  |  .  (2-78) 

2m 

In  obtaining  (2-75)  the  approximation  m  -*■  0  since  m  «  E 

has  been  used.   Equation  (2-75)  is  correct  for  scattering 

from  a  point  proton  with  no  anomalous  magnetic  moment. 

But  the  proton  is  known  to  have  an  anomalous  moment  so 

the  structure  of  Mf.  must  be  changed  to  reflect  this  fact. 

In  addition,  the  possibility  of  nucleon  structure  should 

be  taken  into  account  in  M„.  (the  electron  remains  a  point), 

f  l 

Figure  *1  is  redrawn  to  schematically  illustrate  the  change 
in  Mfi. 
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Figure  5 

The  bubble  at  the  photon-nucleon  vertex  represents  the  non- 
locality  of  the  photon-nucleon  interaction  and  the  anomalous 
moment  of  the  proton.   The  matrix  element  is  now 


M 


fi 


J 


P  q2 


where 


JR -u(pf)oy  U(P.) 


(2-79) 


and 


0y  =  F1(q2)Yy  +  F2(q2)cyvq. 


(2-80) 


The  quantities  F,  and  F„  are  undetermined  real  functions  of 

q  ,  while  a  =  ■*   (y  ,y  }.   Foldy  [25]  has  shown  that  while 

(2-80)  is  not  unique  it  has  the  proper  general  form. 
With  the  above  changes  (2-75)  becomes 


da 

do. 


R 


da 
dft 


x  %   x 


'Mott 


o    *2  2 

4m  c 
P 


4  {^i+vv2  ^  i+4Ft\ 


(2-81) 


where  the  final  bracket  is  the  invariant  matrix  element 


M 


fi 


|  ,  and  where  k   =  1.79  is  the  protons  anomalous  (Pauli) 


52 


moment.   Equation  (2-8l)  is  known  as  the  Rosenbluth  cross 
section.   The  factor  F, (q  )  is  found  to  represent  the 
proton's  charge  distribution  and  its  normal  magnetic  moment, 
while  Fp(q  )  represents  the  anomalous  part  of  the  magnetic 
moment.   F,  and  F?  are  called  the  Dirac  and  Paul!  form 
factors  respectively.   Yennie  [26],  and  later  Hand,  Miller 
and  Wilson  [27]  pointed  out  that  F,  and  Fp  should  not  be 
considered  to  be  electric  or  magnetic  form  factors.   More 
properly  certain  linear  combinations  of  F,  and  Fp  defined 
to  be 


.2  2 

E     1    i,  2  2   p  2 
4m  c 


q2  +  0 


—  1 


(2-82) 


and 


GM  "  Fl  +  W 


2 
q  -*■  0 


—  1  +  <   =  y 
P    P 


(2-83) 


are  termed  electric  and  magnetic  form  factors.   The  major 

advantage  of  (2-82)  and  (2-83)  is  that  in  the  Rosenbluth 

2 
cross  section  (2-8l)  the  cross  term  (F,  +  k  Fp)   is  removed, 


da 
dfi 


da 
dft 


'Mott 


o   ^2  2 
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E   ,.  2  2 
4m  c 


=>/* 


^2  2 


*2  2 

_    *  q 

j,  2  2 
4m  c 


4m  c 


2  tan2  |  G2 


(2-84) 


In  both  the  cross  section  formulas  the  form  of  the 

2        2      2  0- 
quantity  in  braces  is  A(q  )  +  B(q  )  tan   p-%   It  can  be 

shown  that  this  is  the  most  general  form  for  the  one  photon 

exchange  process,  regardless  of  the  structure  of  the 
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electron  and  the  proton.   The  interpretation  of  A  and  B  in 

terms  of  G„  and  G   or  F,  and  F~  depends  of  course  on  the 
E      ml      2   ^ 

structure  assumed. 

2 .   Nucleon  Form  Factors 

In  the  text  above  no  real  attempt  has  been  made  to 
coherently  explain  the  form  factor  concept.   So  far  they 
have  been  characterized  mainly  as  correction  factors  that 
when  multiplied  by  the  scattering  cross  section  for  point 
(structureless)  particles  yield  the  cross  section  for 
scattering  from  physical  particles.   Thus,  as  introduced, 
these  dimensionless  quantities  are  meaningless  unless  some 
predictions  can  be  made  from  a  basic  theory  of  matter  about 
the  structure  and  behavior  of  the  form  factors. 

While  the  simple  idea  of  the  form  factor  representing 
some  physical  extension  in  space  can  be  illustrated  by  non- 
relativistic  examples,  a  more  comprehensive  discussion  of 
their  meaning  needs  the  help  of  the  so-called  dispersion 
relations.   Both  viewpoints  are  investigated  below. 

In  Section  II,  part  A2,  the  scattering  amplitude  f(Q-) 
was  found  in  a  partial  wave  solution  of  the  Schroedinger 
equation  for  the  scattering  of  neutrons  from  protons.   Of 
interest  here  is  the  solution  of  the  Schroedinger  equation 
for  the  scattering  of  electrons  from  two  potentials  V(r), 
one  with  a  point  charge  source,  the  other  with  a  smeared 
out  charge  distribution. 

With  all  assumptions  as  in  Section  II  part  A2,  the 
Schroedinger  equation  for  scattering  is  again 


5^ 


(V2  +  k2Wr)  =  ~V(rH(r).  (2-85) 

This  differential  Schroedinger  equation  is  solved  formally 

for  ^  making  (2-85)  into  an  integral  equation.   The  inverse 

2     2 
of  the  linear  differential  operator  (V   +  k  )  is  an  integrals 

hence 

¥»(r)  =  /G(r,r')  2g  V(r f )*(r ' )dr'  (2-86) 

J  ftd 

where  r*  is  a  dummy  variable  of  integration  and  G(r,r')  is 
a  function  characteristic  of  the  form  of  the  differential 
operator.   G(r,r')  is  known  as  a  Green  Function.   The 
general  solution  to  (2-85)  is  the  sum  of  the  solution  to 
the  homogenous  problem  (in  this  case  it  would  be  the 
unscattered  plane  wave  e     since  —5-  V(r)  =  0)  and  some 
particular  solution  to  the  inhomogenous  problem  which  is 
given  by  (2-86).   Hence 


♦ 


(r)  =  elk'r  +  /G(r,r')  ~   V(r ' )*Kr ' )dr ' .      (2-87) 


This  is  called  the  integral  equation  for  scattering,  and 
the  solution  is  found  by  iteration.   The  first  iterated  form 
is  called  the  Born  or  plane-wave  approximation 

*  =  e1^  +  %  /G(r,r<)  e1^  dr'.  (2-88) 

fi  J 

Comparing  (2-88)  and  (2-16)  yields  the  scattering  ampli- 
tide* 


In  order  that  the  comparison  be  made  the  correct  form 
of  G(r,r')  must  be  determined. 
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f(e-)  =  ^/elq-r  V(r')dr*  (2-89) 

ft  J 


-*■        -> 


where  q  =  k.  -  kf.   For  a  spherically  symmetric  potential 
V,  (2-89)  is 

00 

f(G-)  =  ~pT  /  r'sln  qr'  V(r')  dr'.  (2-90) 

2 
For  a  point  charge  of  Ze,  the  potential  V  is  Ze  /r,  thus 

fpt-(&)  =  Ti4  <2-9D 

^         *h  q 

while  for  a  distributed  charge  of  Ze  of  the  form 

VD  =   5|_   (l-e-ar);    a^C, 

fDW  =  2S|£  _^  =  f       (&)  _^  .  (2.92) 

hqq+a  ^  q+a 

Hence  the  structure  of  the  distributed  charge  shows  up  as  a 
factor  times  the  point  charge  scattering  amplitude.  It  can 
be  shown  that  the  form  factor  here  is  the  Fourier  transform 

of  the  charge  distribution.   Figure  6  illustrates  the 

2 
behavior  of  (2-91)  and  (2-92)  as  a  function  of  q  .   The 

2 
slope  of  f(Q-)  at  q   =  0  is  proportional  to  the  radius  of 

the  structure  under  consideration.   Before  continuing  it 

must  be  stressed  that  all  that  has  been  said  so  far  of  form 

factors  refers  to  the  non-relativistic  case.   In  a  rela- 

tivistic  example  difficulties  are  immediately  apparent.   A 

charge  distribution,  spherical  in  one  reference  frame  is 

elongated  in  another  frame  so  that  the  meaning  of  the  radius 
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f(e-) 


slope  =  0=  structure 
radius  =  0 


slope  i-    0=  finite 

structure 
radius 


Figure  6 

of  a  charge  distribution  is  not  clear.   No  longer  can  a  form 
factor  be  thought  of  properly  as  the  Fourier  transform  of 
the  distribution.   One  has  to  be  concerned  with  the  form 
factors  themselves  and  not  their  three-dimensional  Fourier 
transforms  in  a  particular  reference  frame. 

The  Meson  theory  of  nuclear  forces  mentioned  briefly  in 
Sections  I  part  D  and  II  part  A3  states  that  virtual  pions 
are  constantly  being  emitted  and  reabsorbed  by  a  nucleon. 
Further,  two  nucleons  are  stuck  together  through  an  exchange 
of  pions.   Now,  since  nucleons  are  surrounded  by  a  cloud  of 
pions,  the  virtual  photon  emitted  by  an  electron  scattered 
from  the  nucleon  might  have  in  fact ,  interacted  with  the 
nucleon  meson  cloud.   So  the  simple  electromagnetic  inter- 
action of  a  photon  and  a  nucleon  may  be  quite  a  bit  more 
complex  and  depend  on  the  dynamics  of  the  strong  interaction 
of  the  pion  and  the  nucleon.   Dispersion  relations  then 
attempt  to  understand  form  factor  behavior  as  a  function  of 
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2 

q   in  terms  of  the  strong  interaction  dynamics.   What 

follows  is  a  simple  minded  attempt  to  make  plausible  the 
form  of  the  equations  put  forth  for  the  form  factor. 

The  basic  idea  of  dispersion  relations  (spectral  repre- 
sentations, or  S-matrix  theory  of  strong  interactions)  is 
that  a  quantum  mechanical  amplitude  for  a  physical  process 
is  the  boundary  value  of  an  analytic  function  of  a  complex 
variable.   For  electron  scattering  F(q  )  would  be  the 
boundary  value  along  the  real  line  of  the  analytic  function 
F(z)  where  z  -*  q   +  ie.   A  dispersion  relation  is  an  inte- 
gral equation  that  relates  a  dispersive  process  to  an 
absorptive  process.   The  best  known  example  is  that  of  the 
dispersion  of  light  in  a  dielectric,  where  the  complex  index 
of  refraction  n(w)  =  n  (w)  +  i  -^—   a(w)  is  expressed  as  an 
integral  over  all  frequencies  involving  the  absorption 
coefficient  a(w) ; 


00  oo 

2 n  f     w'lm  n(w)  ,  ,  c  /   a(w')dwT 

=  -  P  /   — ~ ^y*-  dw  '  =  -  /     0  ~ 

tt   /     1 2     2  tt  /     r  2     2 


( 2— 9 3 ) 
r  7T   /     t2     2         77  /     »2     2 


(P  represents  the  Cauchy  principal  value  of  the  integral) 

A  classical  optical  example  of  (2-93)  is  Sellmiers  equation 
for  the  anomalous  dispersion  of  light, 

N    A  A2 
n2  _  i  =  \  *  (2-94) 

i-1  X     ~Xl 

where  X  is  the  incident  light  wave  length,  A.  the  i 
resonant  wavelength,  A.  a  constant,  and  N  the  number  of 
natural  frequencies  of  the  medium. 
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In   electron   scattering 


00 

_1    [  ^ 

T  Jr.        „    ' 


ImP(q    )  0 

F(q^)    =    i       /      — -p 2    dq  C2^95) 

7T 


is  a  direct  analog  of  (2-93).   Now  three  mesons,  the  p,  4> , 
and  the  to-meson  have  been  found  experimentally  to  have  the 
correct  properties  so  that  they  may  couple  to  the  nucleon 
and  to  the  virtual  photon.   They  correspond  to  the  N  natural 
frequencies  above.   Hence  it  is  plausible  to  write  keeping 
(2-94)  in  mind 

F(q2)  =  S     —^-k  j  =  p,  n.  (2-96) 

i-i  i-4 

m. 

i 

In  terms  of  the  physical  form  factors  G,  (2-96)  might  be 

written 

P             gwp             g(f>p             gpp 
GE    (q    )    =    2   +   2    +   2    +    X  "  gcop  "  S*p  "  Spp 

p  x_a_     !_a^     !_a_ 

m  mf  m 

and 

G      (q2)    =      ^P      +   .  y^P      +      ypP      +l,i|       -y       -u       (2-97) 
UM    Kq    '  2  2  2         X      ywp      y<J>p      ypp  ^   y'; 

■"-2  2  2 

m  m,  rn 

w  9  p 

where  the  g's  and  u's  are  coupling  constants  to  be  fitted  to 

experimental  data.   The  twelve  coupling  constants  in  (2-97) 

are  not  all  independent  and  their  number  can  be  significantly 

reduced.   The  best  fit  available  for  the  proton  was  used 

here  to  predict  GE   for  use  in  (2-84).   De  Vries  [28]  adjusted 

P 
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the  free  parameters  appearing  in  (2-97)  so  that  a  minimum 

2 
is  obtained  in  a  x   fit  to  experimental  data.   The  following 

values  were  determined:   g   =  2.628,  g,   =  1.853* 

2         -2     2         —2 
g   =  1.191,  where  m   =  15.6  F   ,  m,   =  26.6  F    are 

experimentally  confirmed  resonances  of  fairly  narrow  width. 

Since  the  p  resonance  appears  as  a  rather  broad  peak  (750  ± 

100)  MeV  and  since  some  observers  [29]  have  shown  that  the 

effective  position  of  the  p  mass  can  be  at  about  600  MeV, 

the  mass  of  the  p  resonance  has  been  taken  as  a  free 

2  -2 

parameter,  determined  to  be  m   =  8.463  F   .   Before  con- 
tinuing it  should  be  noted  that  both  (2-9*0  and  (2-97)  are 
approximations  of  (2-93)  and  (2-95)  respectively  in  the 
sense  that  the  resonances  considered  are  taken  to  be  delta 
functions  in  the  integrands  of  the  more  accurate  form  factor 
expression  (2-93)  and  (2-95) 

It  is  interesting  to  note  that  the  theory  outlined  above 
actually  preceded  the  discovery  of  the  mesons  that  it 
requires.   Dispersion  relations  as  applied  to  the  strong 
interaction  was  developed  in  1957  by  Chew  [30]  et  al  on  the 
basis  of  a  two-pion  exchange  approximation.   This  theory 
was  not  able  to  explain  form  factor  behavior,  and  it  later 
showed  (1959)  that  a  three  pion  resonance  must  exist  in 
order  that  theory  and  experiment  agree. 

A  multiplicity  of  form  factor  notations  exist  today. 
Coupled  with  their  often  misleading  names  and  differently 
defined  static  limits,  they  present  a  baffling  array. 


60 


Commonly,  the  starting  place  is  with  the  two  functions 

-*-2  -*-2 

F, (q  )  and  F?(q  ),  called  the  Dirac  and  Pauli  form  factors 

(really  there  are  four  form  factors,  two  for  the  proton, 
two  for  the  neutron) .   They  describe  the  structure  of  the 
normal  Dirac  coupling  and  the  anomalous  Pauli  coupling 
between  the  nucleon  and  the  electromagnetic  field.   F,  and 
Fp  are  not  directly  related  to  the  electric  and  magnetic 
moment  distributions  of  the  nucleon,  but  in  older  papers 
they  are  often  erroneously  called  electric  and  magnetic 
form  factors.   Still  another  general  name  for  F-,  and  F? 
are  vertex  form  factors  referring  obviously  to  their  use 
in  the  current  operator  J  .   The  Dirac  and  Pauli  form  fac- 
tors have  the  following  static  limits: 

P1  (0)  =1  F1  (0)  =  0 

P  n 

F2  (0)  =  1.79  =  k     F2  (0)  =  -1.91  =  Kn     (2-98) 
p  p      n 

Hand,  Miller  and  Wilson  [27]  popularized  the  use  of  the 
physical  form  factors  G.   They  were  introduced  to  simplify 
the  Rosenbluth  cross  section  by  removing  the  cross  product 
terms  in  F-,  and  F?  and  are  in  no  way  more  fundamental  than 
the  vertex  form  factors.   They  are  defined 

VP1  +  "KF2  "  FCH   ' 

Gm  E  Fx  +  KF2  E  (l+k)PMAG  (2-99) 

■*■?    2 

where  n  =  -q  /km    ,  m  the  neutron  or  proton  mass,  k  the 

anomalous  proton  or  neutron  moment.   F  .  and  F  „   were  used 
^  ch      mag 
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briefly  before  the  G  factors  became  popular.   The  form 

J 

factors  G„  and  G   do  merit  the  names,  electric  and  magnetic 

Em  '  to 

form  factor,  and  they  have  static  limits 

GE  (0)  =1         GE  (0)  =  0 
p  n 

Gm  (0)  =  2.793      Gm  (0)  =  -1.91  (2-100) 

P  n 

It  is  often  found  that  still  another  type  of  form  factor 
is  convenient  to  work  with.   In  isotopic  spin  notation 

where  t  =  1  for  the  proton  and  x  =  -1  for  the  neutron,  the 

z  ^  z  ' 

charge  and  magnetic  form  factors  can  be  written  succinctly 
for  both  the  neutron  and  proton  as 

GE    "    l(1  +  Tz}    GE     +¥1~\)    GE 

P  n 

G      =   1(1  +       )    GM     +1(1  -       )    G 
m        2  zy      M        2  zy      m 

p  n 

or   as , 

GE   =    2(GE     +GE    }    +    2  rz(GE     _GE    } 
p  n  p  n 

G      =    1(G        +.G      )    +   1        (G       -  G      ) 
m        2      m  m  ?   *  z     m  m 

p  n  p  n 

Now  define   the    isoscalar   form   factors    as, 

GE    ~    2(GE      +   GE    } 2 1*  2 

p  n  q      ■>   0 

GS    =   1(G         -   G      ) = -l(u     +y    )  (2-101) 

m        2m  m  2        n      2      P        ny 

p  n  q      -*■    0      ^      ^ 

and  the  isovector  form  factors  as, 
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GE  ~    2(GE   "  GE  }     2     "  2 

p      n     q  -*•  0 

GI  '=    X   "  Gm  > 2TT*2-35  (2-102) 

p     n     q  ■*  0 

These  (iso)  form  factors  are  convenient  for  use  in  theoreti- 
cal work  because  the  electromagnetic  current  J   may  be  split 

s     V 
into  isoscalar  and  isovector  parts,  i.e.  J   =  J    +  J 

H    '      y    y     y 

The  intermediate  states  in  the  photon-nucleon  vertex  can 
be  characterized  by  this  isoscalar-isovector  idea  in  that 

the  two  pion  resonances  with  (T=0,  J=l)  w  and  <j>  manifest 

s       s 
themselves  in  the  isoscalar  form  factors  G^   and  G   , 

E       m  ' 

whereas  the  third  resonance  (T=l,  J=l)p   contributes  to  the 
isovector  form  factors  G^   and  G,ff  .   It  is  interesting  to 

Cj  1v1 

note  that  only  the  isoscalar  form  factor  2GF   =  (G„  +  GF  ) 

p    n 
appears  in  the  deuteron  cross  section,  which  might  be  expec- 
ted since  the  deuteron  current  is  an  isoscalar. 

There  is  an  important  relationship  between  the  physical 

form  factors  that  was  discovered  in  1961  at  Orsay  by 

2 
Lehmann,  Taylor,  and  Wilson  [31] .   In  low  q   measurements 

2         -2 
(actually  the  work  was  at  q   =  2.98  F   )  of  the  e-p  cross 

section  the  following  was  noted: 

n     Gm 
E  -   — 2-  (2-103) 

P   y 
M  p 


_2 
Later  data  verified  this  up  to  approximately  100  F   ,  and 

(2-103)  came  to  be  known  and  called  the  "scaling  rule." 

There  is  some  suggestion  of  (2-103)  when  the  vertex  form 
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factors  are  expressed  in  terms  of  the  physical  form  factors, 
for  example  inverting  (2-99)  yields 

Gt-,  —  T)G  G   —  G-r-, 

n     Em     t-,     m   E 

Fi  =   i-n  •  <F2  =  -TTTT- 

The  above  have  poles  at  n=l  unless  (2-103)  is  true.   Drickey 
and  Hand  have  experimentally  verified  (2-103)  to  include 
the  deuteron  with 
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A  further  relationship  of  this  type  has  been  noted. 


(2-104) 
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Figure  7  shows  that  the  scaling  rules  can  be  used  as  an 
established  empirical  fact  (to  at  least  20  F  )  without 
introducing  appreciable  error. 
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3.   Elastic  Electron  Deuteron  Scattering 

Following  the  same  general  procedure  as  Rosenbluth 
did  with  the  proton,  Jankus  (1956)  derived  [3]  the  cross 
section  for  elastic  electron-deuteron  (e-d)  scattering  from 
a  point  deuteron  in  the  first  Born  approximation.   His 
results  differ  qualitatively  from  the  proton  cross  section 
in  one  important  way,  i.e.,  there  are  three  form  factors 
characterizing  the  deuteron  while  two  suffice  for  the  proton 
This  is  in  agreement  with  the  work  of  Glaser  and  Jaksic  [32] 
Their  study  showed  that  the  cross  section  for  scattering  a 
relativistic  electron  from  a  potential  with  spin  J,  contains 
2J  +  1  form  factors.   The  three  deuteron  form  factors  are 

found  to  correspond  to  the  charge,  quadrupole,  and  magnetic 

2 

moments  in  the  static  limit  q   -*■  0 .   The  general  form  of  the 

cross  section  is  as  it  should  be  for  the  one  photon  exchange 
process,  A(q  )  +  B(q  )  tan   p-. 

Specifically  Jankus  found  that 


££  =  [  d£\     (i         2E    •  2  ± 


L-l 


where 


^2,-^2^    8   2„2,->2x 
Vq  }  +  9  n  Vq  } 


+  |n  |  1  +2(l-n)  tan2  |j  D2(q2) 


oo 


u2(r)  +  w2(r)]jo(i) 


dr 


(2-106) 


(2-107) 


is  the  charge  structure  factor, 
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D  (^  =  6/2   2 
uQKq    }  2d 


0 

I 


w  (r) 


u(r)w(r)-rL-     j2(x)dr 


•5 


M^Q   (2-108) 


is  the  quadrupole  moment  structure  factor. 

CO  oo 

D^(q2)  =  f  [u2(r)  -  |w2(r)]jo(x)dr  +   ~  f   [u(r)w(r) 


+  ^Zllj?(T)d 

•  2   J 


r —  1  -  — P 

r     l    2rD 


(2-109a) 


accounts  for  the  contribution  of  the  intrinsic  magnetic 
moments  of  the  proton  and  the  neutron  to  the  scattering 
process  and 


DM(q  )  =  f  /  w  (r) 


oo 

/  w2(r)[jQ(T)  +  j2(x)]dr 
0 


1  P 

2  D 


(2-109b) 


is  the  magnetic  contribution  to  the  scattering  process 
arising  from  the  convection  of  charge  in  the  deuteron. 
Together  (2-109a3b)  are 


DM<52)  "  | 


?dU  +  dS] 
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M 


M 


(2-110) 


Al 


so  given  are  the  static  limits  (q.  -»-0).   The  j   and  jp  are 


spherical  Bessel  functions,  t  =  |q|r/2,  0  is  the  quadrupole 
moment  and  Pn  the  percentage  D-state  (see  page  127). 

While  the  currents  in  e-p  scattering  were  easily  found, 
the  deuteron  current  presents  somewhat  more  of  a  problem. 
The  virtual  photon  is  expected  to  interact  with  one  nucleon 
only.   But  the  interaction  amplitude  of  an  electron  with  a 
bound  nucleon  is  unknown.   The  natural  thing  to  do  is  to 
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assume  that  since  the  deuteron  is  a  loosely  bound  system, 
the  photon  interacts  with  a  nearly  free  nucleon,  while  the 
other  nucleon  participates  just  as  a  spectator.   This  is 
called  the  impulse  approximation  and  is  described  graphic- 
ally by  Figure  8. 


Figure  8 

Here  the  bubbles  represent  the  photon  -  free  nucleon  vertex. 
Distortion  of  the  free  electromagnetic  form  factor  (vertex) 
due  to  binding,  i.e.,  the  difference  between  bound  and  free 
form  factors,  is  thought  to  be  less  than  5%  [331.   Now  that 
the  impulse  approximation  has  freed  the  nucleons,  the  nucleon 
participating  with  the  photon  in  the  interaction  is  described 
by  the  free  field  Dirac  Equation. 

Since  the  impulse  approximation  involves  the  inter- 
action of  a  photon  with  an  individual  "free"  nucleon,  the 
current  is  again  of  the  form  UO  U,  where  0   is  given  by 
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Op  =  Plp(q  )YM  +  F2p(q  )aM  q^ 

oy  =  p  (q2)yy  +  F  (q2)aPV 
n    In  M   '      2n  M     Mv 


(2-111) 


The  deuteron  current  is  then 


J 


v 


-uoyu  +  voyv 

P      n 


(2-112) 


Adler  [^51  writes  the  relevant  part  of  the  S-matrix  form  in 
the  Impulse  approximation  as 


o    ~  |S,J2  -  | 


e   Ty  |  2 


fi 


M_  I  -  =  M  ,  ~  J 


fi 


y   2  "d1  * 

q 


(2-113) 


Just  as  in  the  proton  case,  | M  . |   is  evaluated  and  found 


to  be 


Mfi|2  =  Gc  +  FGQ  +  K 


1+  2(l-n)  tan2  | 


(2-111) 


+2  /fl„2 


where  n  =  a-q  /^M,,  and  G  ,  GQ  and  GM  are  the  charge,  quad- 
rupole  and  magnetic  form  factors  given  by  (including  their 
static  limits ) . 


Gc  -  (ge   +  ge  >  Dc — r-r1 

p      n       q  =  0 


(2-115) 


and 


and 


GQ  "  (GE   +  GE  >  DQ" 
p     n 


wdQ 


G   =  (G    +  GE  )D*j  +  (G„  +GM  )2DjJ 
p     n         p    n 


(2-116) 


M 


"**  M   yd' 


(2-117) 
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where  the  D's  are  given  in  equations  (107)  through  (110). 
The  steps  from  a  point  deuteron  to  a  finite  size  deuteron 
is  thus  carried  out  by  the  introduction  of  the  free  nucleon 
form  factors.   This  step  is  justified  by  the  impulse 

approximation . 

2  2 

The  contributions  of  G   (charge),  G~    (quadrupole)  and 

2 
GM  (magnetic)  to 

A(Q2)  =  G2(q2)  +  |  n2G2(q2)  +  |nG2(q2)       (2-118) 

are  shown  in  Figure  9. 

h .      Neutron  Charge  Form  Factor  from  e-d  Scattering 

Equation  (2-117)  can  be  simplified  considerably  by 
scaling,  i.e., 

d  d 

and 

GM  GM  n 

_ £  =  _ R  =  GE 

p  n 

For  the    case   of   the    deuteron 

<>M      -    «3E      +    GE    )    »d   DC  (2"119) 

d  p  n 

so   that   with    1   -   n   ->   1, 

„2         ln  „       ,2    r^2/n    L2  i2         i]         t2    ,       2    0-v       8      2n2l 

Gd   =    (GE      +   GE    }         DC(1+3nyd      +    3   ™d      tan      2)  +  9nDQJ 

P  n  L  ~>  _> 

Gd   "    (GE      +   GE    )2    Pd2  (2"120) 

P  n 
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Contributions  to  A(q  ) 


Total 


Figure  9 

The  functions  G   (charge),  —  n  GQ  (quadrupole)  and  =■  nG^ 
(magnetic).   The  smooth  curve  is  the  sum  of  these  three 
terms  (taken  from  Ref.  ^6). 
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where  F,  is  called  the  deuteron  structure  factor  and  is  the 
d 


same  as  the  result  of  Jankus  for  a  point  (nucleons)  in  the 
deuteron.   Prom  (2-120) 

(2-121) 


G,        /    "E 

±Bp'   l+-n 
G^    *d  I    G, 


E 


The  aim  of  the  experiment  is  the  evaluation  of  the 

charge  form  factor  of  the  neutron.   The  magnetic  contribu- 

2       2 
tion  to  G,  and  G   must  therefore  be  eliminated.   In  the  case 
d      p 

of  the  deuteron  the  magnetic  contribution  is  given  by 
(2-120)  as 


c  =  |  ny'2  (l  +2  tan2  £) 


(2-122) 


so  that 
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GEXP 
1  +  C 


with   G 


EXP 


exp 
mott 


d 


(2-123) 


where  the  quadrupole  term  has  been  neglected. 

In  the  case  of  the  proton  the  magnetic  contribution  is 
found  by  considering 


GE  y  GM 

G2   =   — £+nT-£(l+2    tan2   |) 
p       i_n  l-n  2' 


(2-124) 


Scaling  assumes  GM  =  uGF  where  u  =  2.79  so, 


G' 


GE 
=  1^-  {l-^2  (1+  2  tan2  f)}. 


(2-125) 


The  term 


C  =  ny2  (1  +  2  tan2  |) 


(2-126) 
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is  then  the  magnetic  correction  term.   Thus 

gp  "  T^  (1  +  c) 


so  that 


where 


2    Gp<l-1> 

Ge  =  -^ (2-127) 

p    1  +  C 


G2  =  _exp_  < 


P   o 


mott 
Finally,  for  experimental   comparisons 

GE.       /   GE 

GE     Fd  I    GE 
P       \     P 

or 

GE       GE 

p-21  =  ~  pr^  -  1.  (2-128) 

GE     Pd  GE 
P         P 

The  value  of  GF   from  (2-128)  is  very  sensitive  to  the 

n 

uncertainty  in  P,  and  GF  /G„  .   For  example  a  1%   error  in 

n   p 

F  ,  translates  into  about  a  10$  error  in  G„  /G„  . 
d  E   E 

n  o   P 
The  only  reliable  data  in  the  range  of  q^  considered 

here  have  been  presented  by  Drickey  and  Hand  [1].   Their 
results  were  in  complete  disagreement  with  the  neutron- 
electron  interaction  slope.   Their  data  will  be  presented 
after  reanalysis  of  the  theory  in  part  B5 . 
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5.   Review  of  Relatlvistic  Deuteron  Theories 

Essentially,  all  the  material  above  is  based  upon 

the  solutions  to  a  nonrelativistic  Schroedinger  equation  for 

the  two-nucleon  problem.  The  results  produced  using  the 

nonrelativistic  S  and  D  state  wave  functions  are,  at  times, 

surprisingly  good.   Still,  there  are  problems  as  the  nagging 

difficulty  in  the  determination  of  the  percentage  D  state, 

and  the  theoretical  determination  of  the  deuteron  magnetic 

moment,  not  to  mention  the  conflict  that  is  the  point  of 

this  work,  i.e.,  GF   s  0  while  dG,,  /d(q2)   P° 

n  En       Q   =  0- 

Hence  sufficient  reason  exists  for  improving  the  theory, 

as  well  as  the  experiment. 

The  Jankus  theory  on  which  the  cross  section  (2-106) 
was  derived  is  itself  one  of  the  major  limitations  in  the 
development  of  better  answers  for  the  deuteron.   Gross 
states  that  the  theory  gives  no  indication  of  what  has  been 
left  out,  or  how  to  calculate  corrections,  which  might 
include  the  effects  of  of f-the-mass-shell  contributions 
arising  from  the  fact  that  the  nucleons  are  bound  and  not 
free.   Intermediate  states,  i.e.,  pion  exchange  currents, 
could  cause  further  corrections.   Thus  to  isolate  possible 
corrections  a  fully  relatlvistic  theory  should  be  developed, 
one  that  would  reduce  to  Jankus  theory  at  low  q2.   Correc- 
tions then  could  be  localized  and  more  carefully  defined. 
This  is  a  difficult  order  to  fill,  and  to  date  no  entirely 
satisfactory  solution  has  been  proposed,  though  many  authors 
have  tried.   The  most  interesting  attempts  are  briefly  dis- 
cussed below. 
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Van  [3^1]  has  calculated,  using  deuteron  wave  functions 
derived  from  a  Bethe-Salpeter  equation,  relativistlc  cor- 
rections to  the  deuteron  form  factor.   These  wave  functions, 
which  are  relativistic  generalizations  of  Hulthen  wave 
functions  (no  hard  core)  tend  to  make  the  functional  D  of 
Equation  (2-106)  larger  at  high  momentum  transfers  than  do 

models  having  a  hard  core.   The  essential  part  of  Van's 

ch 
results  is  a  comparison  between  his  predictions  of  Dp  (F  . 

in  Van's  notation)  calculated  with  the  (relativistic) 

Hulthen  wave  functions  and  the  corresponding  D~(F  ,  ,  Van) 

of  the  nonrelativistic  case.   The  comparison  showed 

F  ,  (q  )  >  F  ,  in  a  range  of  q   from  0  to  9  F~  ,  which 

results  in  F,  being  larger  than  in  the  relativistic  case. 

Recalling  (2-128)  quickly  confirms  the  fact  that  the  end 

result  is  to  make  GF  /GF   even  smaller.   The  increase  in 

ch       NR        n    P 
F  ,  from  F  ,  is  attributed  to  the  quadrupole  contribution 

2       -2 
(+3.5%  at  q   =  1  F   ),  and  to  final  state  wave  function 

2       -2 
Lorentz  contraction  (-1.5$  at  q   =  1  F   ). 

Adler  [^51  has  obtained  what  amounts  to  first  order 
relativistic  corrections  to  the  nonrelativistic  impulse 
approximation  deuteron  form  factors.   His  theory  is  intrin- 
sically nonrelativistic  and  differs  from  that  of  Section 
II  part  B2  only  by  the  use  of  Paul!  spinors  (for  convenience) 
instead  of  Dirac  spinors,  and  by  retaining  powers  of 

deuteron  moments  such  that  the  final  cross  sections  are 

2   2 
correct  to  order  q  /m  .   The  extra  terms  now  found  in  the 


V* 


deuteron  current  are  incorporated  into  additive  corrections 


to  Gp  and  GQ .   Specifically 
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AGn=2MJ   w(u'r-u)   ~  0 dr  +  — 

'0  qV         5/2. 


r    2ji 

/   w  — 


jn(qr/2)-4j3(qr/2) 
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and  where 


GRS  E  GR   +  GR   »   GR  =  Fl(cl2)  +  2icF2(q2). 
P      n 


At  q2  =  1.0  P  2  where  Gc  =  0.55,  AG^ 


_5 

x  10    and  where 


R  -k 

G0  =  0.03,  GQ  =  2  x  10   ,  hence  the  corrections  are  rather 

too  small  to  effect  F,  greatly  but  go  in  the  right  direction 

Further,  these  corrections  to  F,  do  not  alter  the  value  of 
5  d 

y,  since  it  depends  only  on  P~,  and  the  value  of  Q,  is 
changed  by  a  bare  1% .      In  summary  then,  it  is  doubtful  that 
these  corrections  are  very  meaningful  in  a  theory  using 
nonrelativistic  wave  functions. 

Adler  glossed  over  at  least  one  sticky  problem  in  the 
above  work.   He  admits  the  need  for  a  "fudg'e  factor"  in 
resolving  the  kinematical  difficulty  in  the  lab  frame  of  an 
energy  transfer  (recoil)  to  only  one  of  the  nucleons  by 
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the  virtual  photon.   Gourdin  [351  has  made  what  amounts  to 

the  same  calculation  but  in  the  Breit  frame  where  the 

2    -*2 
momentum  transfer  is  entirely  spacelike,  i.e.,  q   =  -q  , 

2 
q~  =  0.   His  difficulty  was  in  the  transformation  of  the 

nonrelativistic  deuteron  wave  functions  into  the  Breit 
frame.   His  solution  to  this  was  to  assume  that  the  wave 
function  were  unchanged.   Gourdin ' s  results  were  approxi- 
mately the  inverse  of  Adler's,  i.e.,  of  the  same  magnitude 
but  opposite  in  sign.   Unlike  Adler,  Gourdin  did  derive  a 

very  small  correction  to  the  magnetic  form  factor  but  it  is 

2       -2 
so  small  that  it  is  not  observed  for  q   <_  6  F   .   The  value 

then  of  these  works  lies  in  the  fact  that  it  almost  certainly 

seals  the  fate  of  any  straightforward  attempt  to  determine 

(quasi)  relativistic  corrections. 

The  strictly  relativistic  approaches  (a  covariant 

dispersion  relation  representation)  are  more  basic  and  more 

difficult  than  the  quasi-relativistic  work  above.   Success 

with  this  technique  is  limited  to  the  low  momentum  transfer 

region  where  it  is  now  believed  that  the  theory  is  more 

accurate  than  experiment.   The  early  papers  of  Jones  [4]  and 

Cutkosky  [36]  are  interesting  when  applied  to  the  G„ 

n 
problem.   Using  the  triangle  impulse  approximation  diagram 

(page  67)  Cutkosky  calculated  the  simplest  absorptive 

amplitude  which  yields  a  deuteron  structure  factor  of  (all 

spins  ignored) 


t?A/  2s     1 


/TT^     .   -1  /n  C0t  IT 
tan 


/n   cot  jjjj-  /l+n 
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,a=/wM.  (2-129) 


The  relativistic  analog  using  the  tail  approximation  is 


Fd(q2)  =  /  u(r)Jo(qr/2)dr  =  -^-   tan"1  i^J-     (2-130) 


0 


where 


°  "  °Mott  GE2  Fd-  (2-131) 

T       A 
Comparison  of  the  two  F,  and  F,  yield  the  relativistic 

corrections 

T 
T        Fd 
Fd -l4  (2"132) 

and      |q| ■ ^-  (2-133) 

/i+n 

These  imply  that  a  surprisingly  rapid  curvature  is  needed  at 

2 
small  q   for  GE   in  order  that  the  n-e  interaction  be  in 

n  -2 

consonance  with  G^   -Oat0.3-0.6F 

E 
n 

Gross  has  done  the  most  extensive  work  in  the  area  of 
the  purely  relativistic  approach  [2  ,  8,  9,  10,  11].   He 
ambitiously  set  out  to  construct  a  complete  theory  fully 
relativistic  and  fundamental  in  the  sense  that  very  few 
phenomenological  parameters  occur  and  based  on  single- 
variable  unsubtracted  dispersion  relations  and  coupled 
unitarity  equations.   He  had  partial  success  but  numerical 

predictions  of  the  ratio  of  relativistic  to  nonrelativistic 

2 
form  factors  seems  much  too  large  at  high  q  .   In  the  range 

2  -2 

of  low  q   (=0.25-1. OF   )  agreement  is  reasonable,  and 

his  predictions  for  u,  and  Q,  are  within  2%   of  the 
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experimental  value.   In  work  on  relatlvistic  corrections  to 
the  impulse  approximation  in  elastic  e-d  scattering  Gross 
realized  that  the  purely  dispersion  technique  tended  to 
mix  up  the  internal  structure  of  the  interacting  particles, 
i.e.,  the  deuteron  wave  function  and  the  nucleon  form  fac- 
tors with  the  deuteron  current.   This  leads  him  to  eventually 
express  his  theory  in  a  formulation  closer  to  that  of 
potential  theory . 

Specifically  he  obtained  a  modified  Bethe-Salpeter  type 

-2 
wave  function  retaining  terms  up  to  the  order  of  N   ,  and 

then  used  the  relatlvistic  nucleon  current  to  derive  the 
deuteron  form  factor.   An  estimate  of  the  size  of  the  cor- 
rection to  the  charge  form  factor  is  obtained  by  taking  the 

2 
slope  of  the  correction  term  at  q   =0.   Gross  found  that 

p 
the  correction  is  proportional  to  -1/Bm  which  is  of 

right  sign  so  as  to  make  the  neutron  charge  form  factor 
more  positive  (towards  agreement  with  thermal  n-e  slope 
data)  , 

The  relativistic  correction  mentioned  above  came  mainly 
from  relativistic  modifications  of  the  deuteron  wave  func- 
tions and  the  nucleon  current.   Adler  had  previously  looked 
into  the  nucleon  current  modifications  but  his  analysis  had 
been  based  on  nonrelativistic  deuteron  wave  functions. 
Gross  states  that  the  wave  function  modifications  have  the 
largest  effect  and  hence  Adler 's  corrections  are  insufficient. 
Van  did  use  relativistic  wave  functions  from  a  Bethe-Salpeter 
equation  in  his  study  of  corrections  to  the  impulse 
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approximation.   His  difficulty  was  in  using  a  separable 
potential   which  was  equivalent  to  using  a  Hulthen  wave 
function  which  has  no  hard  core,  and  thus  gives  values  of 
F,  too  large.   The  trouble  is  that  here  relativistic  cor- 
rections to  e-d  scattering  and  the  model  dependence 
question  have  been  intermixed.   Van  did  not  avoid  this  by 
comparison  of  his  relativistic  corrections  to  corresponding 
nonrelativistic  models  since  the  two  were  not  in  the  same 
reference  frame.   That  is  to  say  the  relativistic  model  (in 
the  Breit  frame  for  example)  when  Lorentz  transformed  into 
the  deuteron  rest  frame  was  not  compared  to  a  corresponding 
nonrelativistic  model,  which  if  it  is  a  good  model  already, 
simulates  all  relativistic  effects  as  seen  in  the  rest 
frame . 

The  principal  contribution  to  the  correction  of  the 
charge  form  factor  of  the  deuteron  comes  from  taking  into 
account  the  recoil  motion  of  the  deuteron,  which  in  turn 
distorts  the  wave  function.   Lorentz  contraction  also  adds 
to  the  wave  function  distortion,  though  not  as  much  as  in 
Van's  case.   Ignoring  the  correction  for  the  quadrupole  and 
magnetic  terms,  the  result  of  Gross'  calculation  for  the 
charge  form  factor  is 

GC  =  (GE   +  GE  }  DC  (2-134) 

P     n 


at 

Van  did  this  in  order  to  improve  on  Jones  and  Cutkosky's 

e-ar 
asymptotic  "tail"  wave  function  approximation  iK  =  . 
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where 


D 


R 


u+&) 


oo 

/<■ 


2  +  w2)  j   dr  +  ADP 


C2-1.35) 


and 


ADC  = 


a 


2  2 

(u  -  uru'  -  utcu  +  7w  -  wrw '  -  wkw)  (j  +j'  )  dr 


(uru '  +  wrw ' )  j   dr. 
°  o 


(2-136) 


In  this  expression  prime  refers  to  differentiation  with 
respect  to  r,  and 

,2 


k  =  r 


d 


dr 


-  a 


(2-137) 


The  arguments  of  the  Bessel  functions  are  of  course  qr/2, 
and  q  is  |q  |.   Now  in  terms  of  Gross's  corrections,  A(q  )  is 

A(Q2)  =  (G   +GE  )2D^2  =  (GE  +GE  )2(P^+AFd)  (2~138) 
p    n  p    n 

where  the  quadrupole  contribution  to  the  charge  form  factor, 

2  2 
8/9n  GQ,  is  small.   Thus,  since  A  is  an  experimental  quan- 
tity and  does  not  change  due  to  corrections,  a  change  is 

forced  in  G„  .   To  see  this  write 
E 
n 


Gc  =  (QE   +  Gf  )Df  =  (GE   +  Gf  )DC 

p     n  p     n 

(no  corrections). 


(2-139) 
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But 

GC  =  CGE   +  GE   )DC  =  (GE   +  [GER  +  AGE  ^(DcR  +  AV  (2-l4°) 
p     n  p      n      n 

(with  relativistic  corrections). 
Equating  the  right  hand  sides  of  (2-139)  and  (2-1*10)  yields 

-AGE  D^R  -  (GE   +  GgR)  ADC  -  Gg  ADQ 
n         p     n  p 

or 

AD         AF 

AGE   s  ~GE   IW  *  "GE   TT  (2-l4l) 

n      p  Dc       p  M 

where  F,  is  the  deuteron's  charge  structure  function.   Gross 

wrote  (2-1^1)  leaving  out  the  F  , .   Since  F ,  -£  1  for  very  low 

2 

q  ,  AGE  -   GE  AF,  gives  an  upper  limit  to  the  correction  as 

2     n     p 

q  -*■  0  in  the  n-e  interaction.   Figure  10  summarizes  the 

2 
corrections.   In  the  range  of  q   considered  here,  Gross' 

2 

correction  AGE   is  proportional  to  -  -^-p  and  practically 

n  8M  4.     „„u      t  a^   ■ 

model  independent.  p 

Using  this  correction  to  the  charge  form  factor  of  the 

neutron  plus  Feshbach-Lomon  wave  functions  for  the  deuteron, 

Gross  reanalyzed  the  G„   -  0  data  of  Drickey  and  Hand. 

n 
Figure  11  shows  a  comparison  of  results  with  Partovi  (i.e. 

HJ)  and  Feshbach-Lomon  wave  functions.   Gross'  correction 
coupled  with  FL  wave  functions  seems  to  remove  the  discrep- 
ancy between  the  electron  scattering  and  n-e  interaction 
results . 
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The  relativistic  correction  is  labeled  by  AG„   ,  the 

En 

difference  between  the  Feshbach-Lomon  model  and  the  Partovi 

model  is  labeled  GlFL)_  gX?)  and  the  meson  exchange  contri- 
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A  recent  attempt  by  Chilton  and  Urhane  [37]  to  explain 

GE   -  0  for  q2  -  0.3  -  1.0  F~2  with  G^  (0)  /  0,  was  made 

n  n 

by  fitting  the  nucleon  form  factors  with  rational  functions 

2    2 
of  the  variables  1  +  q  /m   .   They  found  G^   is  never  far 

^    it       °  h 

n      2 

from  zero  and  bends  over  rapidly  with  increasing  q   (see 

Figure  10).   Nevertheless  with  Chilton  and  Urhane's  fit, 

_2 
G„   does  not  get  back  to  zero  till  beyond  12  P 

n 
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III.   APPARATUS  AND  EXPERIMENTAL  PROCEDURES 

A.   THE  NPS  LINAC 

The  Naval  Postgraduate  School  linear  accelerator  (LINAC) 
is  a  low  intensity  (20  yamp  maximum  average  current),  100  MeV, 
three  section,  30   foot  long  electron  linear  accelerator. 
While  its  construction  and  general  operation  have  been 
described  elsewhere  [38],  several  important  modifications 
have  been  accomplished  and  will  be  mentioned  briefly. 

A  supplementary  cooling  water  system  with  a  steam  heat 
exchange  was  added  for  the  accelerating  sections.   This  was 
done  so  that  the  water  temperature  could  be  precisely  con- 
trolled (±1°C)  throughout  a  run  lasting  for  an  extended 
length  of  time.   In  the  present  experiment,  runs  of  up  to 
38  hours  were  made,  and  since  the  machine's  energy  stability 
is  strongly  dependent  on  the  accelerator's  changes  in 
temperature,  this  is  a  critical  point. 

Until  recently  the  target  chamber  vacuum  had  been 
isolated  from  the  accelerator  sections  by  a  three  mil 
aluminum  window  32  inches  upstream  from  the  target.   This 
window  caused  considerable  beam  spread  and  by  the  time  the 
target  was  illuminated  by  electrons,  the  beam  spot  was 
approximately  3/16  -  1/k    inch  in  diameter.   Now,  with  the 
addition  of  a  fore-diffusion  vacuum  pump  combination  in  the 
target  area,  the  entire  accelerator,  from  electron  gun 
(source)  to  target,  is  open.   In  fact  with  the  spectrometer 
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coupled  to  the  target  chamber,  electrons  pass  only  through 
the  target  Itself  and  a  five  mil  aluminum  exit  window  from 
the  spectrometer  before  striking  the  detectors.   Upstream 
secondary  emission  monitors  (SEM)  have  been  moved  to  loca- 
tions downstream  of  the  target.   The  accelerator  -  target 

chamber  -  spectrometer  combination  is  maintained  at  a 

_5 
pressure  of  10   mm  of  Hg.   The  removal  of  the  SEM  and 

window  has  significantly  reduced  the  beam  spot  size,  and 

it  is  now  estimated  from  beam  appearance  on  a  zinc-sulfide 

(ZnS)  screen  that  2/3  of  the  electrons  are  confined  to  a 

spot  1.0  mm  in  diameter. 

Past  nuclear  structure  experiments  [391  have  experienced 

severe  problems  in  background  radiation  levels.   Major 

shielding  modifications  were  made  to  correct  this  situation. 

More  paraffin  and  borated  paraffin  in  the  form  of  a  floor 

to  ceiling  wall  was  added  to  the  existing  shielding  in  the 

beam  deflection  area.   In  the  target  area,  a  platform  loaded 

with  paraffin  four  feet  above  the  floor  and  extending  30 

inches  out  from  the  wall  was  suspended  from  the  overhead. 

Its  purpose  was  to  shield  the  counter  house  from  line-of- 

sight  radiation  from  the  beam  dump.   Wax  was  placed  in  the 

spaces  between  the  I-beams  in  the  overhead  in  order  to 

reduce  reflection  of  neutrons  into  the  counting  system. 

The  primary  shielding  effort  was  on  the  counting  house 

itself.   Concrete  blocks  were  replaced  by  boxes  of  paraffin 

and  much  lead  was  replaced  by  special  bricks  composed  of 

5%   boron,  Q0%    lead  by  weight  and  polyethylene.   A  cross 
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section  of  the  counter  house  is  shown  in  Figure  12.   This 
shielding  arrangement  reduced  the  background  to  quite  an 
acceptable  ratio  of  2000  to  1. 

B.   THE  TEN-CHANNEL  COUNTING  SYSTEM 

1 .   The  Ten-Channel  System 

In  general  there  are  two  main  categories  of  detec- 
tors:  slit-defined  momentum  acceptance,  and  counter- 
geometry  defined  acceptance  detectors.   A  single  counter, 
slit  defined  detector  was  replaced  by  a  ten  counter 
geometry  defined  system.   The  ten  counters,  which  span  a 
portion  of  the  focal  plane  of  the  spectrometer  and  define 
adjacent  momentum  intervals,  are  known  as  a  "ladder."   The 
'advantages  of  this  arrangement  over  the  single  channel 
system  is  obvious;  for  one  spectrometer  momentum  setting, 
ten  pieces  of  information  instead  of  one  are  gathered. 

The  ladder  is  composed  to  ten  7/16    inch  high,  1/16  inch 
thick  NE102  plastic  scintillators  placed  approximately  1/2 
inch  behind  the  focal  plane  of  the  spectrometer  as  shown 
in  Figure  12.    A  large  one-piece  scintillator  used  in 
coincidence  with  the  front  counters  and  called  the  backing 
counter,  is  placed  3  inches  behind  the  10  singles  counters. 
The  choice  of  7/16  inch  high  (in  the  dispersion  direction) 
counters  produce  a  momentum  acceptance  Ap/p  of  about  0.3$ 


*The  ten  scintillators  (counters  or  channels)  are  num- 
bered in  ascending  order  1  to  10,  where  channel  1  represents 
the  lowest  momentum  accepted. 
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per  channel  or  about  3%  for  the  entire  ladder.  An  electron 
scattered  from  the  target  and  subject  to  the  spectrometer's 
magnetic  field  B,  follows  a  path  of  radius  r,  given  by 

.  '     r=£=B^-  (3-D 

Thus  higher  momentum  electrons  are  seen  by  the  higher 
numbered  channels.   Further  the  momentum  resolution  is 
greatest  for  the  higher  channels.   An  electron  traversing 
the  spectrometer  and  passing  through  a  particular  scintil- 
lator produces  a  pulse  that  is  collected  by  an  Amperex 
XP  1110  photomultiplier  tube,  amplified,  transmitted, 
discriminated,  fanned  out,  and  finally  registered  on  one 
of  ten  scalers.   A  block  diagram  of  the  counting  system  is 

at 

shown  in  Figure  13.    With  the  exception  of  the  scalers, 
which  allow  a  counting  speed  of  20  MHz,  the  electronics  is 
capable  of  100  MHz  counting  rates. 

2 .   Resolution  and  Efficiency  of  the  Ten-Channel  System 
In  order  that  the  system  be  useful,  two  quantities 
must  be  known:   The  relative  counting  efficiency  of  each 
channel  and  the  momentum  separation  of  adjacent  channels. 
The  momentum  separation  (resolution)  of  the  channels  was 
obtained  by  taking  elastic  scattering  peaks  from  a  thin 
carbon  foil  and  detecting  the  electrons  scattered  as  a 


The  Chronetics  and  EG  and  G  counting  equipment  was 
loaned  to  NPS  by  the  High  Energy  Physics  Laboratory,  Stanford 
University . 
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function  of  spectrometer  momentum  in  each  channel.   That 
is,  the  10  counters  were  treated  as  if  they  were  single 
counters  while  the  sharp  elastic  peak  was  moved  through 
them  by  changing  the  momentum  setting  of  the  spectrometer. 
A  plot  of  counts  in  the  scalers  vs .  spectrometer  momentum 
yields  10  elastic  scattering  peaks  (Figure  14).   The  separ- 
ation of  the  peaks  as  a  function  of  electron  momentum  p  is 
a  measure  of  the  channel's  resolution.   The  10  elastic  peaks 
produced  are  not  identical  in  height  or  area,  but  since 
they  all  "see"  the  same  elastic  carbon  peak,  the  differences 
are  ascribed  to  the  differing  relative  efficiencies  in  detec- 
tion and  recording  of  the  electrons. 

Now  specifically,  the  resolution  of  the  10  channels  was 

determmined  by  taking  21  -  ten  peak  sets  at  incident  elec- 

MeV       MeV 

tron  momentums  from  27  to  9^  and  with  scattering 

c         c 

angles  varying  between  ^5°  and  105°.   These  experimental 
data  were  fit  to  a  quadratic  expression 

p  =  p  (a  +  3N  +  YN2)  (3-2) 

s 

where  N  is  •the  channel  number,  p   is  the  momentum  setting  of 

the  spectrometer,  and  p  is  the  mean  momentum  of  channel  N. 

A  best  fit  produced  values  for  a,  $  and  Y.   If  p   is  defined 

as  the  momentum  setting  for  which  a  particle  of  momentum  p 

entering  the  spectrometer  passes  through  the  center  of 

channel  6,  then  (3-2)  can  be  rewritten  as 

=  p  [1  +  A(N-6)  +B(N-6)2] 


PN(PS)    Fs 
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A  and  B  again  are  constants  independent  of  p.   The  channel 
momentum  acceptance  is  then 


r~^ /P   =  A  +  2B(N-6).  (3-H) 

d(N-6)   /  s 

Typical  values  of  A  and  B  are  3.0  x  10~3  and  -1.3  x  10" 
respectively  and  it  is  clearly  seen  that  the  higher  channels 
(i.e.,  N  >  5)  have  a  higher  resolution  than  the  lower 
channels . 

So  that  the  energy  spectrum  of  the  scattered  electrons 
can  be  determined  properly  the  relative  efficiencies  need 
to  be  measured.   The  efficiencies  can  vary  for  several 
reasons,  some  of  which  are:   Differences  in  detection 
thresholds  and  physical  properties  of  the  plastic  scintil- 
lators; spectrometer  transmission  characteristics;  and 
solid  angle.   The  determination  of  the  relative  efficiencies 
was  based  on  the  method  of  Crannel  and  Suelzle  [^0]  ,  and 
essentially  involved  measuring  the  shape  of  an  inelastic 
continuum  of  scattered  electrons.   The  shape  was  found  by 
observing  the  spectrum  with  all  ten  channels  at  several 
different  momentum  settings  assuming  initially  that  each 
channel  has  efficiency  one.   Then  the  counts  from  individual 
channels  are  compared  to  the  calculated  shape  and  corrected 
channel  efficiencies  are  obtained.   Next,  a  more  accurate 
shape  is  predicted  on  the  basis  of  the  corrected  efficiencies 
and  the  process  is  repeated.   Three  or  four  iterations  will 
normally  produce  constant  efficiencies. 
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Typically,  five  different  momentum  settings  over  a 
total  momentum  range  of  6%    (twice  the  10  channel  span)  with 
at  least  1000  counts  per  channel  per  setting  were  taken. 
It  is  important  that  the  efficiency  data  be  taken  in  the 
energy  range  in  which  the  experiment  is  to  be  performed. 
Hence  if  an  elastic  deuteron  peak  is  located  at  50  MeV,  then 
the  machine  energy  should  be  adjusted  so  that  a  smooth  part 
of  some  inelastic  continuum  lies  at  that  point.   In  carbon 
12,  the  radiation  tail  of  the  elastic  peak  is  sufficiently 
smooth  at  30  MeV  below  the  peak,  thus  for  the  above  situa- 
tion a  machine  energy  of  80  MeV  would  be  set ,  and  the  five 
momentum  (energy)  settings  could  be  taken  at  52,  51 >  50, 
49,  and  48  MeV  for  example.   Carbon  is  obviously  not  suited 
for  the  NPS  machine  because  then  a  practical  limit  of 
65  MeV  would  be  imposed  on  the  incident  energy  for  the 
proton  and  deuteron  peaks  (i.e.,  65  +  30  MeV  =  95  MeV 
needed  for  efficiency  and  any  more  would  exceed  machine 
capability).   Aluminum  was  substituted;  its  inelastic  tail 
contains  peaks  only  at  2.2  and  3.0  MeV.   Arbitrarily  a  value 
of  6  MeV  down  from  the  aluminum-elastic  peak  was  chosen  as 
the  experimental  area  from  which  to  obtain  the  relative 
counting  efficiencies.   Normally  it  is  best  to  exactly  dup- 
licate in  the  efficiency  measurement  the  conditions  of  the 
particular  experiment,  but  this  was  found  to  be  impossible 
in  some  cases.   On  the  high  energy  runs,  94  MeV  (location  of 
elastic  peaks  under  consideration),  94  +  6  =  100  MeV  was  an 
impractical  energy  to  obtain  and  maintain.   Thus  the 
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efficiency  data  were  taken  at  a  somewhat  lower  energy.   Also 
aluminum  not  CH„  or  CD?  was  always  used  for  efficiencies 
since  it  was  desirable  not  to  take  the  risk  of  burning  or 
destroying  the  experimental  targets  during  the  long  high 
beam  strength  used  in  efficiency  runs. 

To  calculate  the  relative  efficiency  of  channel  k  it 
is  necessary  to  predict  the  total  number  of  counts  channel 
k  receives  in  the  L  spectrometer  momentum  settings.   Thus 

(3-5) 


s=l         s=l  p, 


Ph(k,ps) 


Ek(p,Ps)T(ps)dp 


s=l  p£(k,pg) 


where  p  is  the  particle  momentum,  T(p  )dp  is  the  number  of 

s 

electrons  with  momentum  between  p  and  p  +  dp  which  enter 
the  spectrometer,  E  (p,p  )  is  the  efficiency  of  channel  k 
for  particles  of  momentum  p,  and  p,(k,p  )  and  p«(k,ps) 
the  highest  and  lowest  momentum  accepted  by  channel  k. 
Assume  that 


T(ps)  =  ^Rnp 


n+1 

s 


(3-6) 


n=0 


where  R   can  be  determined  (needs  to  be  done  only  once)  by 


n 


10 


a  least  squares  fit  to  the  points  T  (p  )  =   £   C,  (p  )  as 

c   s    k=  n   ^-   s 

a  function  of  p  .   Hence  the  efficiency  of  channel  k  is 


"TV 


given  by  the  ratio 

L 


Ek  = 


Ls=l 


L 


sw/i  svr1 


s=l  n=0 


(3-7) 
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which  is  the  total  number  of  counts  observed  vs.  the  pre- 
dicted total  number  of  counts. 

C.   THE  TARGETS 

Solid  polyethylene  (CHp)   and  deuterized  polyethylene 
(CD  )   were  chosen  as  the  sources  of  protons  and  deuterons. 
The  immediate  difficulty  is  that  these  plastic  targets  would 
be  destroyed  by  heat  (melted)  in  the  concentrated  electron 
beam.   Hence  the  CH„  and  CDp  targets  were  made  to  move  in 
the  horizontal  direction  so  that  the  beam  sampled  a  width 
of  about  one  inch  over  the  target  face.   The  carbon  in  CH? 
and  CDp  has  much  greater  atomic  weight  (12  umu)  than  does 
either  isotope  of  hydrogen  (H  =  1.00783  umu  and  D= 2.01410 
umu)  and  hence  suffers  much  less  recoil  when  struck  by  an 
electron.   Thus  both  the  elastic  proton  and  deuteron  peaks 
ride  well  down  on  the  inelastic  tail  of  their  respective 
carbon  peaks.   In  order  to  isolate  the  proton  and  deuteron 
peaks  it  is  necessary  to  subtract  out  the  carbon  content. 
A  carbon  target  is  needed  then  for  standardization.   Hence 
arrayed  down  the  target  ladder  were  movable  CDp  and  CH? 
targets,  a  stationary  carbon  standardization  target,  alum- 
inum target  for  efficiencies  and  a  ZnS  screen  for  machine 
tuning  and  beam  observation. 

The  CDp  target  was  3  cm  x  10  cm  x  0.10122  cm  and 
specified  to  be  99.9%  enriched.   The  thickness  of  the  CH2 
and  C  targets  were  governed  by  the  desire  to  make  the  number 
of  C  atoms  in  each  of  the  three  targets  as  nearly  the  same 
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as  possible.   With  the  material  available  a  CH„  target 
0,10^3^  ±  0.00070  cm  thick  and  a  0.09398  ±  0.00061  cm  thick 
carbon  target  were  constructed.   The  carbon  target  was 
produced  by  laminating  7  thin  (5  mils)  graphite  foils.   The 
densities  of  the  samples  were  determined  to  be:   p(CDp)  = 
1.081  ±  .006  gm  cm"3,  p(CH2)  =  0.925  ±  0.005  gm  cm"3,  and 
p(C)  =  1.028  ±  0.005  gm  cm"3. 

In  the  course  of  a  run  the  electron  beam  illuminated 
the  CHp  and  CDp  targets  for  as  long  as  8  hours  apiece. 
After  these  long  runs  considerable  browning  of  the  material 
was  evident,  and  the  possibility  of  damage  to  the  target 
because  of  a  change  in  the  composition  of  the  material  was 
considered.   Subsequent  runs  on  the  undarkened  area  versus 
the  darkened  area  gave  no  indication  of  any  decomposition. 
Further  the  darkened  areas  tended  to  fade  in  time  so  no 
permanent  damage  could  be  supposed.   The  browning  of  the 
target  is  assumed  to  be  an  atomic  electron  dislocation 
phenomena. 

D.   COLLECTION  OP  DATA 

After  deciding  upon  the  momentum  transfer  to  investigate, 
the  appropriate  incident  energy-scattering  angle  combination 
was  chosen.   The  basic  measurement  desired  was  the  ratio  of 
the  electron-proton  to  the  electron-deutron  scattering 
cross  section.   Forward  scattering  angles  were  chosen  where 
possible  in  order  to  emphasize  the  electric  interaction. 
With  this  angular  restriction  the  momentum  transfers  ranged 
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?  — ? 

for  the  proton  and  deuteron  from  about  0.1  <_  q   <  0.5  P 
for  the  available  machine  energies. 

The  machine  was  tuned  for  maximum  beam  and  current 
stability.   The  energy  defining  slits  were  set  for 
AE/E  =  0.5%,    and  the  beam  spot  was  carefully  steered  to  the 
target  center  line.   Collection  of  raw  data  consists  of 
recording  the  number  of  counts  shown  in  scalers  1  through 
10,  the  backing  scaler  counts,  the  time  elapsed  while 
counting,  the  spectrometer  setting  in  MeV,  and  the  number 
of  electrons  passing  through  the  target.   All  the  scaler 
data  are  fed  directly  to  a  teletype  machine  that  both  punches 
an  8  track  paper  tape,  and  prints  out  a  hard  copy.   The 
paper  tape  is  then  interpreted  and  IBM  cards  produced  for 
use  in  computer  data  analysis. 

For  this  experiment  data  were  taken  in  the  following 
way.   The  machine  was  tuned  to  an  energy  approximately  6  MeV 
higher  than  the  energy  at  which  the  elastic  proton  peak 
would  occur.   Efficiency  data  using  a  20  mil  aluminum  target 
were  taken.   The  machine  energy  was  then  lowered  to  the 
value  selected  for  the  proton  run.   Comparison  elastic 
carbon  peaks  in  CH?  and  C  were  taken  to  provide  normaliza- 
tion data  used  in  scaling  the  carbon  for  subtraction.   Then 
the  elastic  proton  peak  from  CHp  and  the  underlying  carbon 
tail  for  subtraction  were  obtained.   Machine  energy  was 
again  lowered  to  the  proper  value  and  the  deuteron  runs  were 
completed  in  a  manner  similar  to  that  for  the  proton.   On 
all  peaks  the  counting  rates  in  the  front  channels  were 
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limited  to  75  counts  per  second,  while  the  backing  channel 
was  limited  to  100  counts  per  sec.   Tests  showed  that  the 
counting  system  (presumably  the  scalers  themselves,  not 
the  counting  electronics)  could  digest  counting  rates  of  up 
to  200  counts  per  second  before  any  appreciable  loss  was 
detected.   Since  the  deuteron  is  twice  as  heavy  as  the 
proton,  less  incident  energy  is  required  to  achieve  the  same 
momentum  transfer  as  that  of  the  proton.   In  all  experiments 
a  peak  was  determined  by  three  spectrometer  settings,  hence 
30  points  define  a  scattering  peak.   The  small  energy  steps 
taken  depended  on  the  machine  energy  and  the  relative  energy 
(momentum)  separation  of  the  various  channels.   The  aim  was 
to  achieve  a  uniform  distribution  of  points  over  the  area 
investigated. 

Where  the  peak  is  located  on  the  ladder  is  important. 
The  end  channels,  1,  2,  9>  and  10  tend  to  have  more  widely 
varying  efficiencies  than  the  center  channels  and  they 
experience  edge  effects,  so  they  are  unsuitable  for  the 
peak  determination.   It  was  found  that  placing  the  peak 
initially  between  channels  6  and  7  gave  the  best  results. 
Then  the  top  of  the  peak  was  defined  by  a  single  counter 
reducing  a  potential  problem  of  matching  the  counts  of 
different  channels  at  a  critical  point.   The  peak  height 
was  chosen  on  the  basis  of  time  versus  statistical  accuracy. 
That  is,  a  total  of  4  x  10   counts  are  needed  under  a  peak 
to  insure  0.5%  counting  statistics,  but  the  time  needed  to 
accumulate  that  many  counts  is  often  prohibitive, 
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2 
particularly  at  the  higher  values  of  q  .   It  was  decided 

that  no  more  than  three  hours  per  spectrometer  setting 
would  be  allotted  for  counting. 

In  order  to  extend  the  range  of  momentum  transfers, 
another  set  of  data  was  taken  by  the  author  at  the  Mark  III 
accelerator  at  Stanford  University.   Where  possible,  over- 
lapping sets  of  data  were  taken  to  check  for  any  systematic 

errors  in  the  equipment.   The  range  covered  at  Stanford 

2         -2 
was  0.20  <_  q  <_   0.80  P   .   No  systematic  deviations  between 

the  two  sets  of  measurements  were  found  and  the  data 

reduction  was  performed  independent  of  the  source. 

The  experimental  equipment  at  Stanford  is  basically  the 

same  as  that  at  the  Naval  Postgraduate  School  with  the 

exception  that  a  100  channel  ladder  is  used  and  that  part 

of  the  data  reduction  is  performed  on  an  IBM  7700  on-line 

computer. 
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IV.   DATA  REDUCTION 

A.  ANALYSIS  OP  RAW  DATA 

Raw  data  are  accumulated  in  the  form  of  counts  in  ten 
channels  for  a  single  spectrometer  setting,  and  are  analyzed 
or  "unfolded"  by  matching  the  counts  in  a  particular  channel 
to  the  appropriate  energy.   The  assignment  of  energy  values 
to  channels  is  determined  by  (3-3).   The  raw  counts  are 
first  corrected  for  background,  then  counting  rate  correc- 
tions are  applied.   Finally  the  efficiency  determined  by 
(3-7)  is  used  to  predict  a  corrected  count  for  the  partic- 
ular channel.   In  practice  all  the  above  calculations  are 
carried  out  by  computer.   Typical  unfolded  energy  spectra 
for  both  proton  and  deuteron  runs  are  shown  in  Figures  15, 
16,  17,  and  18. 

B.  CALCULATION  OF  EXPERIMENTAL  CROSS  SECTIONS  AND  FORM 
FACTORS 

The  experimental  cross  sections  for  e-p  and  e-d  scatter- 
ing were  determined  from 


p,d      N 


sc     „   „  2      -1 


0-c        -   Tt 777  K  K,  cm   ster  (4-1) 

Exp   N.   n.  Afi   s  b 
^     in   t 

where  N   is  the  number  of  scattered  electrons  detected  by 
sc 

the  spectrometer,  N   is  the  number  of  incident  electrons, 

in 

2 
n,  is  the  number  of  target  atoms  per  cm  ,  Afi  is  the  solid 

angle  subtended  by  the  spectrometer  entrance  port,  and 
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finally  K   and  K,  are  the  radiative  corrections  of  Tsai  I4l] 
s       o 

and  Bethe-Ashkin  [H2]  . 

At  the  heart  of  the  calculation  of  (4-1)  is  N   which  is 

sc 

the  integrated  area  under  the  scattering  peak  in  counts-MeV 
per  MeV.   The  scattering  peak  is  produced  by  drawing  a  smooth 
curve  through  the  data  points  in  the  counts  vs.  energy  plot. 

Both  the  elastic  proton  and  deuteron  peaks  ride  on  the 
radiation  tail  of  the  elastic  carbon  peak  (see  Figure  15, 
etc.).   Thus,  in  order  to  obtain  the  relevant  areas  for 
the  cross  section,  the  underlying  carbon  tail  must  be  sub- 
tracted out.   Comparison  of  the  elastic  carbon  peaks  in 
CHp  and  CD?  to  that  of  the  pure  carbon  target  provided  the 
target  normalization  factors.   These  factors  consisted  of 
(1)  "shifts"  of  the  elastic  carbon  peak  from  the  carbon 
target  such  that  its  peak  energy  coincided  with  that  of  the 
carbon  peak  from  CH   or  CD? ;  and  (2)  vertical  "scaling"  of 
the  carbon  tail  through  a  comparison  of  the  areas  of  the 
elastic  carbon  peaks.   The  area  remaining  after  the  sub- 
traction is  N 

sc 

The  number  of  electrons  incident  upon  the  target  is 
measured  by  a  Secondary  Emission  Monitor  (SEM)  downstream 
from  the  target.   The  SEM  is  not  an  absolute  device  to 
measure  the  charge.   In  a  ratio  experiment,  the  only 
requirement  is  that  it  be  stable  during  a  set  of  runs. 
This  stability  has  been  frequently  checked  against  a 
Faraday  cup.   Obviously  the  number  measured  is  wrong  by  at 
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least  the  number  of  electrons  scattered  by  the  target,  but 
this  Is  so  small  compared  to  N.   as  to  be  totally  Insig- 
nificant,  n^  is  a  static  parameter;  it  depends  only  on 
target  composition  and  target  angle  <f> , 
pt  N 


n,  =  

t    A  cos  <J> 


(1-2) 


-3 


where  p  is  the  target  density  in  gm  cm   ,  t  thickness  in  cm, 

N  Avogadro '  s  number,  A  the  atomic  weight,  and  <f>  the  angle 
a 

between  the  beam  line  and  the  downstream  target  normal.   The 

solid  angle  Afi  is  defined  by  1  inch  thick  tungsten  slits  on 

the  spectrometer  entrance  port,  16  inches  from  target  center. 

It  can  be  seen  by  comparison  of  two  experimental  cross 

sections  that  the  only  remaining  quantities  that  enter  a 

ratio  experiment  are  N    and  n,  . 
^  set 

As  mentioned  above,  we  have  used  an  expression  due  to 

Tsai  to  account  for  the  radiative  correction  K  .   This  cor- 

s 

rection  is  a  significant  improvement  over  that  given  by 
Schwinger  [43].   Presented  here  is  Schwinger's  much  simpler 
expression  (4-3)  as  an  illustration  of  the  parameters 
entering  the  correction. 


v  s   r     2a 
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The  Bethe-Ashkin  correction  is  given  by 
6, 


-(l  +  < 


Ale" 
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(4-4) 


107 


Briefly,  K   is  a  multiplicative  correction  to  the  cross 
s 

section  which  accounts  for  radiation  emitted  during  a  large- 
angle  scattering  event  (nuclear  bremsstrahlung) .   K   corrects 
for  the  radiation  that  occurs  before  or  after  the  large-angle 
event  and  is  sometimes  known  as  thick  target  bremsstrahlung. 
For  a  discussion  of  the  terms  in  (4-3)  and  (4-4)  see  Gordon 
[44].   The  E,   used  above  corresponds  to  the  nucleon  recoil 
factor  in  (2-77) . 

The  ratio  of  the  experimental  cross  section  to  the  cross 
section  predicted  for  a  point  nucleon  (Mott  cross  section) 

Equation  (2-76),  yields  the  experimental  form  factors  G 

2  P 

or  G,.   To  these  quantities,  the  magnetic  corrections  as 

outlined  by  equations  (2-122)  and  (2-126)  are  then  applied. 

This  yields  the  square  of  the  electric  form  factor  of 

2 
the  proton,  G„  ,  and  the  sum  of  the  squares  of  the  electric 

P  p       Q    p    O 

and  quadrupole  form  factors  of  the  deuteron,  G^   +  —  r]      Gn  . 

Ed   9     Qd 
Because  of  the  smallness  of  n  in  our  range  of  momentum 

transfers  (n     =  2.2  x  10~3  at  q2  =  0.80  F~2)  the  quadru- 
max  ^  n 

-4 
pole  term  contributes  only  about  one  part  in  10    to  the 

form  factor  and  can  be  safely  neglected. 

The  end  result  of  the  experiment  is  then  the  ratio 

2    2 

G^  /G^   from  which  the  ratio  G„  /G„   is  extracted  with  the 

dp  n   p 

help  of  equation  (2-128). 

In  summary,  the  quantity  evaluated  from  the  data  is 

GL        Nd   x  nP  x  Cd  x  qP  .  .  x  Kd  x  Kd 

d  _   sc t Mott £> b  /j.  c-n 

G2     Np   x  nd  x  Cp  x  ad  , ,  x  Kp  x  Kp 
E      sc     t         Mott     s     b 
P 


10? 


where  C  represents  the  magnetic  correction,  aM  .,  the  Mott 
cross  section,  and  the  indices  p,d  refer  to  proton  or 
deuteron,  respectively. 

C.   EXPERIMENTAL  ERRORS 

Errors  encountered  in  the  measurements  are  of  both 
statistical  and  systematic  nature.   Systematic  errors,  such 
as  the  value  of  the  solid  angle,  counting  rate  correction, 
target  angle,  scattering  angle,  and  SEM  efficiency  cancelled 
out  since  they  were  measured  under  the  same  experimental 
conditions.   The  absolute  error  in  the  energy  of  the  incom- 
ing electrons  is  about  0.2$.   Since  the  Mott  cross  section 
for  the  deuteron  and  the  proton  are  evaluated  at  nearly  the 
same  energy,  this  error  becomes  insignificant.   The  same  is 
true  for  the  errors  in  the  radiative  corrections  if  they  are 
evaluated  over  roughly  the  same  energy  cut  off  AE  (Equation 

4-3). 

The  only  significant  error  of  a  systematic  nature  is 
the  error  in  the  target  thickness  and  density,  pt .   Because 
of  the  thinness  of  the  targets,  it  could  be  established  to 
only  0.5%. 

The  statistical  error  consists  of  two  parts,  the  actual 
counting  error  and  the  error  in  the  efficiency  correction. 
It  is  in  principle  possible  to  keep  these  errors  arbitrarily 
small  by  taking  a  very  high  number  of  counts.   Restrictions 
in  counting  speed  as  posed  by  the  electronic  equipment  and 
time  limitation  in  the  use  of  the  accelerator  pose  a 
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practical  limit.   All  but  two  peaks  contain  40,000  counts, 
which  would  make  their  statistical  error  to  be  0.5$.   The 
required  carbon  subtraction  introduces  an  additional  error, 
depending  on  the  number  of  carbon  counts  under  the  peak  area. 

The  error  in  the  efficiency  corrections  is  quite  ini- 
formly  0.5%  for  the  Naval  Postgraduate  School  runs,  where 
roughly  the  same  number  of  counts  have  been  taken  to 
establish  the  relative  efficiencies.   In  the  Stanford  runs 

the  error  varies  from  0.27%  to  0.65%. 

2    2 
For  each  measurement  of  GF  /Gp  ,  the  error  is  then  the 

d   p 
square  root  of  the  quadratic  sum  of  the  counting,  efficiency, 

and  target  error. 

One  half  of  this  percentage  error  enters  the  ratio 

Gp  /Gp   and,  therefore,  the  calculation  of  GE  /GR  .   This 

dp  n   p 

is  the  error  quoted  in  Table  VII. 
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V.   COMPARISON  TO  THEORY  AND  DISCUSSION 

A.   EXPERIMENTAL  RESULTS 

The  experimental  data  are  presented  in  Table  VI.   The 

error  quoted  consists  of  counting  errors  (including  the 

carbon  subtractions  and  efficiency  corrections)  and  an  error 

of  0.5%   in  the  target  thickness. 

Table  VII  shows  four  sets  of  values  of  Gp  /GF   as 

n   p 
obtained  by:   (a)  the  Feshbach-Lomon  wave  function  with 

relativistic  corrections;  (b)  the  Feshbach-Lomon  wave  func- 
tion without  relativistic  corrections;  (c)  the  Partovi 
wave  function  with  relativistic  corrections,  and;  (d)  the 
Partovi  wave  function  without  relativistic  corrections. 

Values  for  the  relativistic  correction  AGF   and,  as  an 

n 
example,  the  structure  factor  F,  for  the  Feshbach-Lomon 

wave  function  are  also  given.   In  each  of  the  columns  are 

shown  the  slopes  (dGF  /GF  )/dq   to  illustrate  the  rapid 

n   p 
decrease  in  slope  when  going  from  a  relativistically 

corrected  Feshbach-Lomon  wave  function  to  an  uncorrected 

Partovi  wave  function. 

The  b'  fit  of  de  Vries  [28]  for  the  absolute  proton 

charge  form  factor  was  used  to  extract  GF   from  the  GR  /GR 

n  n   p 

ratio.   This  fit  is  in  very  good  agreement  with  the  absolute 

measurements  of  GF   by  Drickey  and  Hand  [1] .   Table  VIII 

P 
shows  the  results  of  the  b'  fit  for  the  two  extreme  cases 

presented  in  Table  VII.   The  slopes  as  determined  by  the 
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the  two  sets  of  data  are  also  given  for  comparison  with 
the  neutron-electron  interaction  slope. 

B.   DISCUSSION 

The  Figures  19  and  20  illustrate  the  results  tabulated 
in  Tables  VII  and  VIII  respectively.   It  can  be  seen  from 
the  slope  of  Figure  20  that  the  Partovi  wave  function  with- 
out relativistic  corrections  does  not  significantly  differ 
from  the  data  of  Drickey  and  Hand  pictured  in  Figure  10(a). 
That  is,  they  both  essentially  agree  that  the  slope  is  very 
small  or  zero.   On  the  other  hand  the  slope  for  the  Feshbach- 
Lomon  wave  function  with  relativistic  corrections  approaches 
closely  the  value  of  the  slope  for  the  neutron-electron 
interaction . 

Barring  unexpected  fluctuations  in  the  proton  form 

factor  G^   in  the  range  of  momentum  transfers  of  interest 
E 
P 
here  it  can  be  concluded  that: 

(a)  The  Feshbach-Lomon  wave  function  together  with  the 

relativistic  corrections  AGF  ,  removes  the  apparent  discrep- 

n  2 

ancy  between  the  neutron-electron  slope  at  q   =0  and  the 

slope  given  from  values  of  G^   (obtained  by  electron  scat- 

n         -2 
tering)  in  the  range  0 .  10  <_  q2  <_  0 .  80  F 

(b)  Even  within  the  relatively  large  errors  propagated 

into  G^  ,  the  Partovi  wave  function  with  relativistic 

E  ' 
n 

corrections  applied,  is  in  disagreement  with  the  neutron- 
electron  interaction  slope. 
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(c)  With  the  present  small  experimental  uncertainties, 
it  is  reasonable  to  state  that  this  experiment  should  be 
added  to  the  group  of  experiments  which  must  be  explained 
by  any  deuteron  theory.   That  is,  a  proper  model  for  the 
deuteron  must  be  found  that  will  predict  the  slope  as 
determined  by  the  neutron-electron  scattering  experiments. 

The  important  qualitative  conclusion  from  these  experi- 
mental results  is  that  the  neutron  has  a  non-zero  charge 
form  factor.   This  implies  that  there  is  a  charge  distribu- 
tion within  the  neutron.   The  neutron  thus  must  have  a 
layered  charge  structure  with  the  outermost  part  of  the 
charge  distribution  being  negative  in  sign. 
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TABLE  VIII 


SLOPE  OF  THE  CHARGE  FORM  FACTOR 

'aq 
'n 


dGE  /dq' 


q2 
(F-2) 

Gi7    a 

E 

P 
de  Vries 

bf    fit 

PFL 
n 

FL 

GE   +AGE 
n           n 

GE 
n 

0,10 

0.9881 

-0.0027 

-0,0020  ±  0,005^ 

-0.0038  ± 

0.0054 

0.20 

0.9765 

+0,0039 

+0.0051 ±  0.0031 

+  0.0017  ± 

0.0031 

0.30 

0.9652 

-0.0006 

+  0.0012  ±  0.0035 

-0.0037  ± 

0.0035 

0.40 

0.9540 

+0.0077 

+  0.0103  ±  0.0036 

+  0.0036  ± 

0.0036 

0.50 

0.9431 

+0.0037 

+  0.0067  ±  0.0074 

-0.0015  ± 

0.0074 

0.60 

0.9321 

+0.0011 

+  0.0047  ±  0.0066 

-0.0042  ± 

0.0066 

0.80 

0.9H5 

+0.0138 

+  0.0184  ±  0.0059 

+  0.0076  ± 

0.0059 

Slope : 

0.0179  ± 0.0036 

0.0026  ± 

0.0036 

Slope  of  n-e  Interaction  =  0 .  0193 ±  0 . 0004 


a:   No  error  has  been  applied  to  Gg  . 
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APPENDIX  A 

THE  NEUTRON  CHARGE  FORM  FACTOR  FROM  THE 
NEUTRON-ELECTRON  INTERACTION 


In  the  n-e  experiments  what  is  measured  is  a  scattering 
amplitude  of  a  neutron  by  a  bound  electron.   In  the  Born 
approximation  for  neutrons  of  very  long  wavelength  this  is 
related  to  the  volume  integral  of  V(r),  the  interaction 
potential   [25] •   The  magnitude  of  the  n-e  interaction  is 

expressed  in  the  form  of  an  equivalent  potential  Vn  over  a 

e2 
sphere  of  radius  rn  =  p,  the  classical  electron  radius, 

u   m  c^ 
e 

so  that  the  volume  integral  of  V„  is  the  experimental  value 
of  the  potential.   The  potential  V"n  is  simply  a  convention 
since  the  classical  electron  radius  plays  no  fundamental 
role  in  the  problem.   Thus, 

jV(r)dr -  Vn   /      dr  =  ?f  rn3  V„    (A-l) 

J  V(r)    =   VQ    =    const.       U  JQ  5       U      U 

Using  Poisson's  equation  and  the  definition,  of  r.m.s.  radius 

2       2      -*■ 
<r  >  =  /r  p(r)dr,  (A-l)  becomes 

V(r)dr  =  -^r—  /r  p(r)dr  =  — - —  <r  >  (A-2) 

Equating  the  right  hand  sides  of  (A-l)  and  (A-2)  yields,  in 
rationalized  MKS  units 

27Te2         <r2>    =    ^-  r  V      (eV)    1.6  x  10~19 


3  X  ^7T£  3    0  0 

o 

2 
where  e      is  the  permitivity  of  free  space  and  in  (Fermi) 
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„    V  (eV)  x  r3  x  8ttg 
<r  >  =  ^-^ F 


e  x  10 


30 


<r2>  =  3.1^  xlO  5  V   (eV)  F2 . 


(A-3) 


From  (2-90)  the  expression  for  the  charge  form  factor  is 
given  by 


2<  2> 
G,-,   =  1  —  -z. +  ... 

E  6 

n 


(A-4) 


so  that  the  slope  at  the  limit  q   +0  gives 


<r2> 


-5.23  xl0~6  V  (eV)  F2      (A-5) 


Krohn  and  Ringo's  [16]  latest  value  is  V  =  -3720  eV,  hence, 


+0.0193  P  . 
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APPENDIX  B 
PROPERTIES  OP  THE  TENSOR  OPERATOR  S 

Several  experimental  facts,  like  the  nonadditivity  of 

the  magnetic  moments  (i.e.,  y   +  y  i-    y,)  and  the  presence 
o  '   p    n    d  v 

of  the  quadrupole  moment  Q,  in  the  deuteron,  led  to  the 
suggestion  that  the  nature  of  the  nuclear  force  is 
partially  noncentral,  that  is,  it  contains  a  partly  "tensor" 
force . 

Require  that  the  tensor  operator  be  symmetric  in 
particles  1  and  2,  and  that  it  be  invariant  under  the  parity 
transformation  (space  reversal).   Let  r  =  x..  -  Xp  be  the 
relative  coordinates  of  particle  1  and  2.   The  particles 
have  spins  given  by  a,  and  a?  where  a's  are  the  Pauli  spin 
matrices  (an  explicit  representation) 


0    1  \        /  0   -i 


These  matrices  obey  Pauli  matrix  algebra, 

[0±i    o.}_   =  26lj     i,j  =  x,y,z  (B-2) 

oxo   =    2ia  (B-3) 

a2  =  a2  =  a2  =  1  (B-4) 

x    y    z 

(a-J)(a-S)  =  (I-S)  +  io-(XxS)  (B-5) 

Application  of  the  parity  operator  P  to  r  and  c  yields, 
Pr  ■>  -r,   Pa  -*■  a. 
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Thus  the  vector  r  can  appear  only  an  even  number  of  times 
in  the  tensor  operator.   Further,  since  the  tensor  operator 
is  symmetric  in  1  and  2,  a,  and  c^  should  appear  in  some 
bilinear  combination.   With  (B-2,3,^,5)  in  mind  it  is  found 
that  the  only  linearly  independent  quantities  which  can  be 
made  from  r,  a,  and  c?  are 

(a^'Op)        and      (a,*r)(a2-r)  (B-6) 

The  traditional  definition  then  for  the  tensor  operator 
S_  p  is 

S12  =  3(a1-r)(c2'r)  -  (c^-Op)  (B-7) 

where  r  =  r/ | r | . 

The  tensor  operator  (B-7)  has  some  interesting  proper- 
ties, namely  <S1?>  over  all  directions  vanishes.   Consider 


=  ihf(°i 


<(a    -r)(o    -r)>    =  ■—-  J(a    -r)(a    -r)    sin  0- d6   d<|>         (B-8) 

1  2  Ht\  J       1  l 

Take  o„    along  z-axis ,  o^  in  the  x-z  plane  making  angle  a 

with  the  z  axis,  and  r  a  vector  making  angles  9-,  and  <J> . 
Then 

a1  =  a1z  (B-9) 

a?   =   <jp(sinax   +  cos  az)  (B-10) 

r  =    sin   0-  cos  <J>  x  +   sin   9-  sin  <j>  y    +   cos   9-  z  (B-ll) 

Putting    (B-9, 10, 11)    in    (B-8)    and   integrating   over  9-  and   <|> 
yields 


122 


<(a1«r) (a2*r)>  =  -  a1a2  cos  a  =  -  (a-^c^) 


(B-12) 


Hence  it  follows  that  <S12>  =  0. 
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APPENDIX  C 

CALCULATION  OP  DEUTERON  MAGNETIC 
AND  QUADRUPOLE  MOMENTS 


The  neutron  and  proton  magnetic  moments  y  ,  and  y  have 

&  n s      p 

been  quite  accurately  determined,  as  has  the  deuteron's 

magnetic   moment    y,.      The    difference    (y      +    y    )    -    u  ,    =    0.0222nm 
o  d  p  n  d 

is  small  but  not  zero  and  much  larger  than  the  experimental 
error.   A  possible  explanation  is  that  the  deuteron  spends 

a  small  part  of  the  time  in  some  state  other  than  the 

i. 
symmetric   S,  state.   The  task  is  to  find  the  correct  mixture 

of  states  that  can  account  for  the  deuteron's  moment. 

Recall  the  angular  wave  functions  for  the  ^S,  and   D. 


states 


3vYoii  =  Y°o*i  <2-^ 


3D    ■   V       1   =J2  Y2    y"1    -f^Y1   y°    +  U-Y°   v1      (2-57) 
ur   7  211      \5     2  xi       \10     2   Xl    \10     2   xl      Kd   ou 

Using  (2-56)  and  (2-57)  construct  the  deuteron  wave  function 
*£SjCr).   Thus 

3v  *oii(?)  ■  ^Voxi  <c^ 


The  nuclear  magnetic  moment  operator  can  be  written 
z  A 

K   -    I   (*zk  +  *pSzk>  +  *n  I       Szk  (C"3) 

k=l  k=z+l 


12^4 


where  JL  is  the  orbital  angular  momentum  of  the  k   proton, 
S,  the  spin  angular  momentum  of  the  k   proton  (or  neutron) 
and  where  g   and  g   are  the  gyromagnetic  ratios  of  the 
proton  and  neutron,  g  =  5.59,  g  =  -3-83.   (The  carot  over 
y   denotes  y„  as  an  operator.)   In  the  c-m  system  and  for 
the  deuteron  (C-3)  becomes 

y      =i&      +gS        +gS  (C-^) 

z        2      z        top    zp        ton    zn 

where  the  orbital  angular  momentum  associated  with  the 
proton  is  half  of  the  relative  orbital  angular  momentum. 
The  magnetic  moment  is  found  by  taking  the  expectation  value 
of  (C-4)  between  the  state  <J>noT  • 


y  >  =  /  < 
z   J 


,-*■    ,  #m      ^  ,m           ^,m      ia    ,    i                         f         . 

*£SJ    yz    *£SJ  B<*ASjlVz"*ASJ>  (C"5) 

The  following  formula  are  useful  in  evaluating  (C-5) 

111  Q      1      1    1 

Szpxl  =  2  Xl  Sznxl  =  2  Xl 

*zp  X?  -  \  A  SznX0  .  -  1  x0         (C-6) 

-1       1    -1  Q      v"1         2   v-1 

S_  Xn   =  -  n  Xn  S     Xn    "  "  ~  X- 


zp  Al      2  Al  zn  Al       2   1 

0  _  1   0  0  _    1   0 

zp  xo  "  2  Xl  zn   !      2  xl 
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Hence 


Oil 


1 
211 


I  (gp  +  gn> 


1 
Oil 


\   (gp  +  ^K   A 


\\- 


gp+gn 


Y2  x 


i-/-L  Y1 

10   2 


0^./      \   ol 

_x1+(g  -gn)  x0J 


sffi 


&n  Y0   1 
2     2  Xl' 


(C-7) 


The  orthonormality  of  the  spherical  harmonics  Y  integrated 
over  the  solid  angle,  and  the  orthonormality  of  the  spin- 
wave  functions  xc   applied  to  (C-7)  yields 


si-  <<J)oiilyzl(f)oii>  =  "^ 


^/ 


u 


(r)dr  =   -E- 


n 


(C-8) 


=    0.8797  nm 


TV  <^;iyJ^)>=  [f  -  \  (gp+gn) 


I.  a 

'211|yz' Y2: 


/ 


w    (r)dr   = 


(C-9) 


J-}  (gp  +  gn}]    "    °-3101 


nm 


In  evaluating  the  radial  integrals  the  assumption  was 
made  that  the  deuteron  was  in  a  definite  state,  either  S  or 
D  but  not  both.   But  the  deuteron  magnetic  moment  is 
0.8573  nm,  that  is,  neither  of  the  two  values  above.   Thus 
consider  a  mixture  of  the  two  states. 

The  total  wave  function  for  a  combination  of  states  is 


*  =  *s  +  ^D  =  ♦oil  +  *211 


(C-10) 
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Normalization  gives 

» 

«{>|<j»  =  Pg  +  PD  =  1  (C-ll) 

co    2  o° 

where  P   =  /   u  (r)dr,  Pn=  /   w  (r)dr  are  the  individual 

S    0  D        0 

probabilities  that  the  deuteron  is  in  an  S  or  D  state.   Thus 

y3s+3D   =   m>+«n>P8+[f-T<«p+Sn>]PD     (C_12) 

Equations  (11)  and  (12)  are  two  simultaneous  equations  in 

the  two  unknowns  P   and  Pn.   Using  the  experimental  value 

s      u 

of  \x      for  u  in  (C-12)  the  solution  for  P   and  PD  are 

PD  =  k%>      Ps  =  96% 

The  quadrupole  moment  Q,  can  be  found  in  the  same  manner 

as  the  magnetic  moment, 

~2    ~2    /      2      s  2  / 

Qd  =  3z   -  rd   =  (3  cos^e-  -  l)r  (C-13) 

In  the  c-m  system,  the  distance  of  the  neutron  and  proton 
from  the  center  of  mass  is  r/2 ,  and  only  the  proton  con- 
tributes to  Q  , .   Hence, 

<Qd>  =  r-/<J>*(3  cos2e-  -  lHr2dr  (C-l4) 

In  terms  of  spherical  harmonics  this  becomes, 


J b2%A  A%V ^\  A%\ 


<Q  >  =  i  JUi 
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3 

The  first  term  represents  the  pure   S..  contribution  and  when 
integrated  over  the  solid  angle  it  is  zero  because  of 
spherical  symmetry.   The  next  terms  in  parentheses  are 
hermitian  conjugates  of  each  other  and  the  total  can  be 
written  as 


dr 


1    fl6rf    f     <?J  *    1   v0  %!     1 

2\  —  JUW7011    Y2   7211 


which   after   angular   integration   gives 

2 


1      f 

— -  /    uwr      dr 

RD  J. 


•50    0 

The  third  term  comes  from  the  JD_.  state  and  gives 


20  J 


2  2, 

w  r   dr 


Hence 


CO  CO 

<Qd>  =  ~J  uwr2  dr  -  ±J    w2r2  dr  (C-15) 

0  0 

To  go  further,  explicit  solutions  for  the  radial  wave  func- 
tions u(r)  and  w(r)  are  needed.   Some  examples  of  Qd  for 
specific  numerical  solutions  u  and  w  are  listed  in  Table  III. 
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APPENDIX  D 
THE  BREIT  OR  BRICKWALL  FRAME 

The  Breit  frame  is  a  Lorentz  frame  in  which  either 

->-    ->  i      ->■    ->- , 

k.,  +  kp  or,  p.  +  Pp  vanish,  where  in  an  elastic  two-body 

-»--»-  ->- » 

collision,  k-,  ,  p.  are  incident  particle  3-momentum  and  k?, 

p'  are  outgoing  particle  3-momentum.   The  somewhat  odd 

notation  stems  from  traditional  techniques  in  elastic 

scattering  theory  and  dispersion  relations,  particularly 

the  Mandelstam  representation.   Consider  Figure  D-l.   The 

convention  is  that  all  four-momentum  are  ingoing.   This  has 

the  advantage  that  one  may  consider  any  two  of  the  four  to 

be  the  incoming  particles  and  the  other  two  as  outgoing 

ones;  the  physical  momentum  of  the  outgoing  particles  is 

then  the  negative  of  the  one  used  in  Figure  D-l. 


Figure  D-l 

The  physical  momentum  of  an  outgoing  particle  is  denoted 
by  a  prime.   For  example,  if  the  particle  corresponding  to 
the  arrow  with  kp  is  outgoing,  it  has  physical  momentum 

k2  "  -k2- 
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Apply  a  Lorentz  transformation  to  the  scattering  event 

such  that  k   +  k~  =  0.   Thus,  four-moment urns  kn  and  k~  have 
1     2  12 

the  form  (n  =  c  =  1) 


k1  =  (w,k)  ,    k2  =  (oo,-k) 


(D-l) 

From  conservation  of  energy  it  follows  that  the  energies  of 
the  p  particle  before  and  after  the  collision,  e,  and  e~ 


must  be  equal.   Hence  ]p-.|  =  |p?|  and 
Px  =  UjP-^  ,    p2  =  (e,p  ) 


(D-2) 


with 


P- 


=  |p. 


=  P 


■V 


2     2 
e   -  m 


Since  scattering  is  elastic  k,  -  k~  =  p?  -  p,  and  it  follows 
k1  -  k2  =  (0,2k)  =  p2  -  px  =  (O^-p^)        (D-3) 

2k  the  three  momentum  transfer.   Equation  (1)  (2)  and  (3) 
yield  the  following  picture.   Both  particles  seem  to  be 
reflected  on  a  hard  wall,  the  k  particle  perpendicularly  and 
the  particle  at  some  angle 


/  Brick  Wall 
/ 

p\  +  p' 


Figure  D-2 
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In  the  c-m  system  the  energies  c    and  w  of  the  p  and  k 
particles  are  not  independent,  but  they  are  in  this  system, 
Therefore  they  are  convenient  to  use  as  variables.   In 
this  system  the  four  momentum  transfer  is  very  simple: 

q   =  (k1  -  kp)   =  (0,2k)   =  -  |2k|   = -q 

i.e.,  it  is  the  three  momentum  transfer. 
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